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ABSTRACT 

We present internal surface-brightness profiles, based on HST/ACS imaging in the 
F606W bandpass, for 131 globular cluster (GC) candidates with luminosities L ^ 10"*- 
3 X 10^ Lq in the giant elliptical galaxy NGC 5128. Several structural models are fit to the 
profile of each cluster and combined with mass-to-light ratios from population-synthesis mod- 
els, to derive a catalogue of fundament al structural and dynamical pa ramete rs parallel in fo rm 
to the catalogues recently produced by i McLaughlin & van der Marell and by lBarmby et alJ for 
GCs and massive young star clusters in Local Group galaxies. As part of this, we provide 
corrected and extended parameter estimates for another 18 clusters in NGC 5128, which we 
observed previously. We show that, like GCs in the Milky Wa y and so me of its satellites, 
the majority of globulars in NGC 5128 are well fit by isotropic 'Wilson models, which have 
intrinsically more distended envelope structures than the standard King lowered isothermal 
spheres. We use our models to predict internal velocity dispersions for every cluster in our 
sample. These predictions agree well in general with the observed dispersions in a small num- 
ber of clusters for which spectroscopic data are available. In a subsequent paper, we use these 
results to investigate scaling relations for GCs in NGC 5128. 
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1 INTRODUCTION 

The spatial structures and internal stellar kinematics of old glob- 
ular clusters (GCs) contain information on both their initial con- 
ditions and their dynamical evolution over a Hubble time. An ef- 
ficient way of extracting this information is to fit detailed models 
to the surface-brightness profiles and (where available) velocity- 
dispersion data of individual clusters, and then to look for possi- 
ble correlations between the physical properties of large numbers 
of GCs. This has been done for most of the 150 globulars in 
the Milky Way, yielding compr ehensive catalogues of their struc- 
tural and dynamical parameters ( Diorgovski 1993; Pryor & Meylan' 
Ll993; Harris 1996; McLaughlin & van derMarel 2005). Much of 
this work has traditionally started from the assum ption that in- 
dividual globulars are well described by the classic I Kind imdt) 
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models of single-mass, isotropic, modified isothermal spheres, 
although recently alternative mode ls have also been employed 
jMcLaughlin & van der Marelll2005h . 

Explorations of numerous scaling relations and inter- 
dependences between the properties of Galactic GCs (e.g., 
iDiorgovski & Mevlanlll994{) have led to the definition of a fun- 
damental plane for globulars that is analogous to bu t physically 
distinct from that for early-type galaxies and bulges jPiorgovskil 
ll995l ; lBurstein et al.|[T997l ; lMcLaughliij |200Q'). Understanding the 
GC fundamental plane in full detail is still not a completely solved 
problem, but important advances have been made in recent years 
as it has become possible to measure the internal properties of GCs 
in many other galaxies. High-resol ution Hubbl e Space Telescope 
(HST) imaging has been used to fit lKing| (ll966^ and other models 
to the surface-brightness profiles of scores of globulars in the Large 
and Small Magellan ic Cloud s and the Fornax dwarf spheroidal 



2003 jIbR iMcLaughlin & van der Marell 



jMackev & Gilmorg 

. M31 ( e.g..lBarmbv Holland. & Huch ra 2002; Barm bv et all 
iioOTI) . M33 jLarsen et alj l2002h . and the giant elliptical galaxy 
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NGC 5128 = C entaurus A ( iHolland, Cote. & HesseJ 1 19991 : 
iHarris et alj l2002h . Internal velocity dispersions and dynami- 
cal mass estimates are also available for smalle r but g r ovying 
numbers of GCs in these systems cDjorgovski et alj 
bubath & Grillmair 19971: iLarsen et al.ll2002l : iMartini & Hoi 
Reikuba et al.,.2007i) . 



1997 



2004 



Here we add to this database with structural measurements 
of 131 GCs in NGC 5128. This galaxy is an attractive target for 
su ch studies in part be cause of its large GC population, estimated 
by [Harris etZI ( l2006l) at Ngc - 1500. It thus contains many ob- 
jects at the high end of the star-cluster mass range (lO'-lO^ A^o), 
which is largely unprobed in the ten-times smaller GC system of 
the Milky Way but where it is increasingly suggested that the clus- 
ter population may encompass a variety of objects including clas- 
sic globulars, the compact nuclei of dwarf elli ptical galaxies, and 
the new class of ultra-compact dwarf galaxies jHilker et al.ll 19991 ; 
[Drinkwater et al. 2000; Hasegan et al. 20o3). In addition, the prox- 
imity of NGC 5128 (D = 3.8 Mpc; see below) makes it possible to 
resolve the core radii as well as just the half-light radii of GCs over 
nearly their full mass range (M > 10* Mq), and thus to fit them 
rigorously with detailed structural models. 

This paper is the third in a serie s of four dea l ing wi th HST ob- 
servations of GCs in NGC 5128. In [Harris et all ilOOj = |PaDerj) 
the Space Telescope Imaging Spectrograph (STIS) and Wide Field 
Planetary Camera 2 (WFPC2) aboard HST were used to measure 
surface-br ightness profiles for 27 very bright GCs in NGC 5128. 
iKind ( [1966) models were fitted to these profiles to derive a struc- 
tural fundamental plane th at could be compare d directly to that of 
the Milky Way globulars. In lHarris et al. l( l2006l = |Paperll we pub- 
lished the first results from a new HST-based survey of GCs in NGC 
5128, using the Advanced Camera for Surveys (ACS) in its Wide 
Field Channel (WFC) to image a total of 131 GC candidate s at a 
resolution of 0'.'05 (linear resolution ^ 0.9 parsec). IPaper III gives 
the full description of the cluster sample along with some rough 
overall characteristics of the ensemble of objects. In the present pa- 
per, we derive surface brightness profiles for all of these clusters 
and fit each of them with a n umber of different struc tur al models . 
The final paper in this series jMcLaughlin et al.]|2007l = |Paper ivl) 
uses these results to examine a number of structural correlations 
for GCs in NGC 5128, which are then compared to the globulars 
in the Milky Way and to var i ous ot her types of massive clusters. In 
related work, iBarmbv et alj ( |2007|) present the results of a similar 
ACS study of GCs in M3 1 and compare the fundamental planes of 
old GCs in that galaxy, the Milky Way, NGC 5128, the Large and 
Small Magellanic Clouds, and the Fornax dwarf spheroidal. 

In the next Section we describe the steps we have ta ken to de- 
rive surface brightness profiles for the GC candidates from lPaper III 
to characterise the point-spread function (PSF) that blurs the very 
central regions (R < 2-3 pc) of these profiles, to transform the 
surface-brightness data from their native HST filter to the standard 
V bandpass, and to estimate metallicities for the clusters from sep- 
arate, ground-based Washington photometry. 

In ^ we apply publicly available population-synthesis mod- 
els to estimate individual V-band mass-to-light ratios for the clus- 
ters, given their metallicities and assuming various (old) ages. Fol- 
lowing this, we summarise the main properties of ea ch o f three 
structural models (those of lKing|ll966llwilsonlll975l and lSersid 
1 19681) that we have convolved with the ACS/WFC PSF and fit to 
every observed surface-brightness profile. 

SectionOgives the results of these fits and uses them to infer a 
wide range of structural and dynamical parameters, including total 
cluster luminosities and masses, effective and core radii and stel- 



lar densities, concentration indices, relaxation times, total binding 
energies, predicted central veloci ty dispersions, and ^-space param - 
eters for the fundamental plane iBender. Burstein. & Fabe j | 19921) . 
We present these in tables that are available in machine-readable 
format either onlin^ or upon request from the first author. Note 
that our measurements of GC luminositie s and intrinsic s izes i n 
particular supersede the recent estimates of Ivan den Berghl ( l2007h . 
who based his numbers on a less detailed a nalysis of some very 
basic cluster characteristics given in jPaper ij 

Injfflwe also address the question of whether the standard 

I Kind ( 1 19661) model specifically gives the best possible fit to GC 
surface-brightness profiles in NGC 5 128. We then extend the range 
of physical parameters calculated for the smaller sample of clus- 
ters previously fitted with iKingt models in .Pa per I. and we provide 
important corrections to some of the more basic parameters (in par- 
ticular, the intrinsic central surface brightnesses) already published 
in that earlier work. The results of this re-analysis are tabulated in 
Appendix IaI 

In ^ we combine our structural modeling and population- 
synthesis mass-to-light ratios to predict line-of-sight velocity dis- 
persions within a series of circular apertures with physical radii 
suited to realistic observational set-up s. We compare these pre- 
dictions with spectro scopic data from iMartini & Hoi ( |2004|) and 
iReikuba et al. for some of our current cluster sample. Fi- 

nally, ^summarises the paper. 

Our modeling analysis in this paper is in all respec t s very 
similar to that undertaken bv lMcLaughlin & van der Marell ( |2005|) 
for a sample of 103 old GCs and 50 young massive clusters 
drawn from the Milky Way, the Large and Small Magellanic 
Clouds, and the Fornax dwarf spheroidal. The catalogues of struc- 
tural and dynamical properties that we produce here for GCs in 
N GC 5128 are likewise very cl ose in form and content to those 
m |m cLaughlin & van der Marell . We have recently completed the 
same type of modeling and produced parallel catalogues for a fur- 
ther 93 G Cs in M31 (Barmbv et al. 20Q7J. As we mentioned above, 
iBarmbv et al.l combine results to co mpare the fundamental planes 
of the old GCs in all six galaxies. In IPaper IvI jMcLaughlin et al.] 
|2007|) we directly compare GC structural correlations only between 
the Milky Way and NGC 5 128, but we also examine how they relate 
to other kinds of massive star clusters. 

In all of what follows, we adopt a distance of 3.8 Mpc to NGC 
5128. This value is representative of recent measurements based 
on the tip of the red-giant branch [(m - M)o = 27.98 ± 0.15], the 
planetary nebulae luminosity function [(m - M)q = 27.97 ± 0.14], 
surface-brightness fluctuations [(m - M)q = 27.78 ± 0.10], Mira 
variables [(m - M)o = 21.96 ± 0.11], and Ceph e ids [(m - M)o = 
27. 67 ± 0.20]; se e'Harris, Harris. & Pooll ( ll999l) . lReiku"bal ( |20()4|) . 
and lFerrarese et al. (2007). The nominal average of these five, rea- 
sonably high-precision distances is (m - M)o = 27.88 ± 0.06, or 
3.76 ±0.11 Mpc. All these methods have undergone rece nt cali- 
bration revisions of various kinds (cf. lFerrarese et al.l2007l) but the 
net results have been to shift the mean up or down by amounts at 
the level of only 0.1 mag. At a distance of 3.8 Mpc, 1 arcsecond is 
subtended by 18.4 parsec. One ACSAVFC pixel (0'.'05) then corre- 
sponds to 0.92 pc. 



See http:// www.astro.keele.ac.uk/~deni/clusters.html | 
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Table 1. NGC 5128 clusters measured independently on two fields " 



Cluster 


Fields 


Cluster 


Fields 


(1) 


(2) 


(1) 


(2) 


AATl 18198 


C018, COW 


C171 


C007, C025 


AAT 120976 


C007, C025 


C173 


C007, C025 


C007 


C007. C025 


C176 


C007, C025 


C018 


C018. CO 19 


F1GC20 


€007, C025 


C025 


C007, C025 


G221 


C007, C025 


C104 


C007, C025 


G293 


C007, C025 


cise* 


C018, CO 19 


PFF021 


€003, C030 


C158 


C018, C019 


WHH22 


€018, C019 



" Boldface in column (2) denotes the field in which the cluster in column 
(1) is closest to the centre of the chip. Model fits to the intensi ty profiles 
from these images are the ones used in the correlation analyses of lPaper ivi 
* Possible star; see Table|2] 

2 DATA 

The GC sample from lPaper ij consists of 62 previously known clus- 
ters in NGC 5128, and 69 newly discovered candidates. All these 
objects fall in 12 target fields imaged in the F6Q6W ("wide V") 
band on the ACS/WFC. We observed 16 clusters twice, since they 
appeared in two overlapping target fields. These are listed in Ta- 
ble[T] We measured two independent surface-brightness profiles for 
each of them, so that in all we have 147 profiles for 131 distinct ob- 
jects. In 94.3l we use these duplications to assess whether variations 
in the PSF over the ACS field of view might have systematically af- 
fected our results. 

Of the 27 GCs observed with STIS or WFPC 2 in lPaperll 9 
were re- observed with the ACS/WFC for lPaperlll and this paper; 
these are C007, C025, C029, C032, C037, C104, C10 5, G221 and 
G293. Our analysis of them here supersedes that in IPaper ll We 
eventually fold the othe r 18 STIS/ WFPC2 clusters into the sample 
for correlations work in lPaper f vl, although with structural param- 
eters updated as discussed in g4.6l and Appendix IaI below. 

We repeatedly convert between luminosities and masses by 
assigning individual K-band mass-to-light ratios to all GCs in our 
total sample. As we describe in more detail below ( g2.3l and 93.11 . 
to do this we first estimate a metallicity for each cluster from 
its (C - Ji) colour in the Washington filter system, using a rela- 
tion calibrated against genuinely old GCs. Then, we input this and 
an assumed old age (normally 13 Gyr) to a standard population- 
synthesis code. Spectroscopy indicates that most of the GCs in 
NGC 5128 are indeed old, but a younger (few Gyr) population can- 
not be ruled out at this stage. If some of the objects in our sample 
are young, then the mass-to-light ratio we assign to them would be 
slightly high, and all physical parameters deriving from it would be 
slightly biased. One way to guard somewhat against this is not to 
include exceedingly blue objects with unrealistically low inferred 
metallicities (see Table |5]l in detailed studies of parameter correla- 
tions and the like. 



2.1 Surface-brightness profiles 

We have used the STSDAS ELLIPSE task to obtain F606W 
surface-brightness profiles for all cluster candidates from lPaper ij . 
As part of the same HST program (GO- 10260), we also obtained 
ACS images for a series of clusters in M31. These data were re- 
duced and mod eled simultaneously with the present sample and 
are discussed in iBarmbv et al] l l2007h . Full details of the surface- 



photometry procedures are given in that paper Here we note that 
we forced the isophote ellpticity in ELLIPSE to be identically 
at all radii. We thus always have circularly symmetric If6 06(R) pro- 
files, which we then fit with spherical structural models. In lPaper 3. 
we showed that the actual ellipticities of the clusters are generally 
quite small, averaging (e) = 0.08 over our whole sample. The as- 
sumption 6 = for the purposes of modeling is therefore not a 
significant limitation. 

The raw output from ELLIPSE is in terms of counts per sec- 
ond per pixel, which we convert to cts/s per square arcsecond by 
multiplying by 400 = (1 px/0'.'05)". Normally, these counts would 
then be transformed to F606iy surface brightnesses, calibrated on 
the VEGAMAG system according to (ACS Handbook) 

y"«06/mag arcsec"^ = 26.398 - 2.5 log(cts/s/n") . (1) 

However, we quickly found that the average, global sky background 
that was automatically subtracted from each ACS image during the 
multi-drizzling in the data reduction pipeline often underestimated 
and sometimes overestimated the local background level around in- 
dividual clusters. The latter case in particular led to the occurrence 
of pixels with unphysical negative counts. We thus had to work im- 
mediately in terms of linear intensity (which we chose to express 
right away as Lq pc"^) rather than going through the usual loga- 
rithmic yU/r606 and then converting to linear quantities later. For the 
solar magnitude we adopt A/q f 506 ~ 

4. 640 and combining this with 

equation l[T) gives 

lF606/Le pc-- - 0.8427 X (cts/s/n") . (2) 

To begin with, we obtained iFsoeW profiles out to R > 10" 
(more than 180 pc) for all clusters. This limit exceeds the expected 
tidal radius for most of them, but such a large field of view en- 
ables us to correct for the inaccurate average sky subtraction in the 
multi-drizzling, by fitting the profiles with PSF-convolved struc- 
tural models that include a constant background term (allowed to 
be negative). We did have to exclude a number of isophotes from 
most of the intensity profiles during this fitting, however. 

First, at very large radii the signal from most clusters is clearly 
swamped by noise, and thus we restricted all fitting to radii R < 150 
WFC pixels, corresponding to R < 7'.'5 - 140 pc. 

Second, the central pixels in some of the brighter clusters were 
saturated. We adopted a saturation limit of 70 cts/s/px and did not 
fit to any isophotal intensities brighter than this. This corresponds 
to a "good" data range of about If606 < 2.36 x 10* Lopc"-, or 
Mf606 > 15.28 mag arcsec"^. 

Third, at intermediate clustercentric radii there are some in- 
dividual isophotes with ELLIPSE intensities that deviate strongly 
from those of immediately neighboring isophotes. To prevent such 
"blips" from skewing the model fits, we first ran the ELLIPSE out- 
put through a boxcar filter to make a smoothed cluster profile, and 
identified points deviating from this by more than twice their own 
(internal) uncertainty. Such points were not included in the error- 
weighted model fitting of ^ 

Fourth, the ELLIPSE estimates of isophotal intensities at clus- 
tercentric radii R < 2px = Qf.'l are all derived from the data in the 
same innermost 13 pixels; but the task nevertheless outputs bright- 
nesses for 15 radii inside 2 px. Clearly not all of these are statisti- 
cally independent. In order to avoid having such correlations (and 
excessive weighting of the central regions of the cluster) bias our 
fits, we decided to include only the ELLIPSE intensities reported 

^ See |http://www.ucolick.org/~cnaw/sun.htinT] 
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Table 2. Irregularities in some NGC 5128 cluster profiles 



Cluster Comments 



AATl 1 1563 Bright star nearby. Fits restricted to i? < 5'.'5 ^ 107 pc. 

AATl 13992 Profile dips at R = I'.'S (30px), near image edge. Fits restricted to R < 1" ^ 20 pc. 

AATl 18198 Diffraction spike from nearby star. Fits restricted toR < 2". 

C104 ON C025 Two nearby stars. Fits restricted to R < 5" ^ 97 pc. 

CI 18 Bright object nearby. Region 0'.'55 <R < 1" masked out of fits. 

CI 34 Bright star nearby. Fits restricted to i? < 3". 

CI 37 Bright object nearby. Region 1'.'2 < i? < 2'.' 2 masked out of fits. 

C154 Cluster C153 aiR^ 2'.'5. Model fits restricted to R < 1" ^ 20 pc. 

C162 Profile dips slightly around R ^ 2'.'5, but fits not biased. No restrictions. 

CI 68 Next to very bright star. Fits restricted to R < 1", which likely misses some cluster. 

Object excluded from sample for correlation analyses in Paper IV . 
C171 Stars at R = 3'.'5 and R = 6'.'6. Fits restricted to R < 6'.'6 ^ 130 pc. 

C174 Star at R = 1'.'6. Fits restricted to R < 1" ^ 20 pc. 

F1GC34 Bright object at R = 2'.'5, and F1GC14 at R = 3'.'3. Fits restricted to R < 1". 

Object excluded from sample for correlation analyses in PagerlV. 
F2GC14 In middle of image artifact (bright star ghost?). Fits restricted to R < 0'.'7 ^ 13.5 pc. 

Object excluded from sample for correlation analyses in PagerlV. 
G170 Profile dips due to image edge nearby. Fits restricted to R < 2'.'1 ^ 40 pc. 

WHH22 Diffraction spike from nearby star at R = 4'.' 8 (star also in sky annulus). 

Fits restricted to R < 3'.'7 ^ 70 pc. 

C145 Very compact, not clearly resolved. Possible star. 

CI 52 Very compact, not clearly resolved. Possible star. 

CI 56 Very compact, not clearly resolved. Possible star. 

CI 77 Very extended; half-light radius R;, =: 6 kpc in a jCmg Q966) model fit. Possible background galaxy. 



Table 3. 147 ^6061^ intensity profiles for 131 GCs in NGC 5128 



Name 


Detector 


Filter 


R 


If606 


uncertainty 


Flag 








[arcsec] 


[Le pc-2] 


[Le pc-2] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


AATl 11563 


WFC 


F606 


0.0260 


1984.157 


19.537 


OK 


AATl 11563 


WFC 


F606 


0.0287 


1963.700 


18.380 


DEP 


AATl 11563 


WFC 


F606 


0.0315 


1939.018 


17.086 


DEP 


AATl 11563 


WFC 


F606 


0.0347 


1907.682 


15.106 


DEP 


AATl 11563 


WFC 


F606 


0.0381 


1875.142 


13.808 


DEP 


AATl 11563 


WFC 


F606 


0.0420 


1836.267 


14.456 


DEP 


AATl 11563 


WFC 


F606 


0.0461 


1792.333 


15.835 


DEP 


AATl 11563 


WFC 


F606 


0.0508 


1741.324 


17.563 


DEP 


AATl 11563 


WFC 


F606 


0.0558 


1680.360 


16.994 


OK 


AATl 11563 


WFC 


F606 


0.0614 


1612.842 


15.518 


DEP 


AATl 11563 


WFC 


F606 


0.0676 


1538.833 


14.016 


DEP 



A machine-readable version of the full Table |3] is available online 
(http://www.astro.keele.ac.uk/~dem/clusters.html I or upon request from the 
first author. Only a short extract from it is shown here, for guidance regarding its 
form and content. Note that the reported f606W-band intensities are calibrated 
on the VEGAMAG scale, but not corrected for extinction. In terms of magnitude, 
^FdOi = 0.308 for the average foreground E{B - V) = 0. 1 1 in the direction of NGC 
5128; see 32.3l for more details. 



for the innermost unsaturated isophotal radius, = i^niin (which was 
always at least 0.5 px), and then for {R-^^ + 0.5 px), 4-1.0 px), 
(^min + 2.0 px), and all R > 2.5 px = a' 125. 

Finally, we looked at every intensity profile individually and 
in a number of cases found irregular features, which we masked out 
by hand. These are summarised in Table[2] 

At the end of Table|2]we also note three GC candidates (C145, 
C152, and C156) that are probably foreground stars rather than 
clusters in NGC 5128, and one (C177) that is more likely a back- 



ground galaxy. We retain these objects in our catalogues of intensity 
profiles and model fits, but only for completeness; none is included 
in any physical analyses. 

Table[3]gives our final, calibrated F606W intensity profiles for 
the 131 objects in our sample (including the duplicate profiles for 
the 16 in Table [TJ. These are not corrected for extinction, which 
we discuss below. Note that only the first few lines of the table are 
reported here; an ascii file containing the full data can be obtained 
from the first author or online. Most of the columns in this table 
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Figure 1. F605W intensity profile, from Table [3] of an object that is not 
clearly resolved. Open squares represent points flagged as "OK" in the ta- 
ble; asterisks are points flagged as "SAT", "DEP", or "BAD". The dashed 
line traces the PSF in equation ^3), after adding a constant background and 
normalising to run through the innermost unsaturated datapoint. 



are self-explanatory (the second and third, which are always WFC 
and F606, are presen t only for compatib ility with the analogous 
table for M31 GCs in lSarmbv et alj|2007l where the detector and 
filter vary from cluster to cluster). The final column gives a flag for 
every point, which can take one of four values: "BAD" if the radius 
is beyond our upper limit of 7'.' 5 or the intensity value is otherwise 
deemed dubious according to the third or final points just above; 
"SAT" if the isophotal intensity above our imposed saturation limit 
of ^ 23, 600 Lq pc"-; "DEP" if the radius is inside < 2 px = 0".' 1 
and the isophotal intensity is dependent on its neighbours (as per 
the fourth point above); or "OK" if none of these apply and the 
point is used when we fit models. 



2.2 Point-spread function 

The ACSAVFC has a scale of a'05 = 0.92 pc per pixel, and thus 
most globular clusters (with typical effective radii Ri, ~ pc) are 
clearly resolved with it. Their apparent core structures, however, 
are still strongly influenced by the point-spread function. Rather 
than attempt to deconvolve the data, we instead fit structural models 
after convolving them with a simple analytic description of the PSF. 
From .Barmbv et al., (2007,) , for the combination of the WFC and 
FmeW filter, 



/p 



f//o = [l + (i?/a'0686)' 



(3) 



which has a full width at half-maximum of FWHM = (X'125, or 
about 2.5 px. Since this PSF formula is radially symmetric and the 
models we fit are intrinsically spherical, the convolved models to 
be fitted to the data are also circularly symmetric. 

Figure [T] illustrates the shape of this PSF and compares it to 
the intensity profile from Table[3]for C156, one of the three GC 
candidates noted in Table[2]as being potential stars. 



2.3 Extinction, transformation to standard V, and cluster 
metallicities 

Once we have fitted models to our clusters' brightness profiles, we 
will want to correct the inferred intensity/magnitude parameters for 
extinction, and to transform the ACS/WFC magnitudes in F606VK 
to standard V for easy comparison with catalogues of other old 
GCs. At the same time, we are interested in predicting observable 
dynamical properties of the clusters (e.g., projected velocity dis- 
persions), which will require some estimate of a mass-to-light ratio 
to apply to our surface-brightness fits. We use population-synthesis 
models to predict V-band MIL ratios, and these require as input an 
assumed (old) age and an estimate of [Fe/H] for every cluster. 

The effective wavelength of the F6Q6W bandpass is 
X^s = 5917 A . Using this in the formula developed by 
ICardelli. Clayton, & MatM ( Il989h implies the extinction relation 
A f606 - 2.8 X E{B — V ), in good agreement with the conclusions 
of lSirianni et alj ( l2005h . The foreground reddening in the direction 
of NGC 5128 is E{B - V) = 0.11 mag, which we adopt for all ob- 
jects in our sample; thus, A^aog = 0.308 mag in all cases, and the 
intensities in Table [3] all n eed to be multiplied by a corrective fac- 
tor 1.328. As discussed in lPaper il the only object for which this 
might seriously be in error is the cluster CI 50, which is projected 
near the central dust lane in NGC 5128. 

Transforming the extinction-corrected intensities to standard 
V is a two-step process, [sirianni et al. I ( 120051) give two transforma- 
tions from F6061V to V magnitude, both including a linear depen- 
dence on de-reddened (V - i?)o colour (see their Table 22). Over 
the colour range 0.4 < (V - R)q < 0.6, appropriate to old glob- 
ulars with [Fe/H] < 0, these two transformations are offset from 
each other by about 0.05 mag. We therefore take their average: 
(V - F606)o = -0.042 + 0.461(V - R)q, with an estimated rms 
scatter of ^ 0.03 mag. Reassuringly, we find that this transformation 
agrees very well with the relation between (V-F606)o and iV -R)q 
predicted for old glob ulars in the population-synthesis models of 
iMarastonI ( fT99ll2005l) . [The VEGAMAG (V - F606)o colours in 
this model were kindly computed for us by C. Maraston.] 

However, we have {C -T\) colours on the Washingto n photo- 
metri c system for most of the clusters in o ur sample dHarris et al.l 
1 19921 ; iHarris. Harris. & GeisleJ |2004|; see IPaper llh . rather than 
(V - R) indices. GeisleJ 1 1996h gives an accurate transformation 
between (V - R)o and (C - T[}q (see his Table 4), and combining 
this with our (V - R]o-(V - F606)o relation yields 



(K- F606)o = -0.019 + 0.126(C- ri)o , 



(4) 



for which we estimate a precision of about ±0.05 mag. We em- 
phasize again that this conversion is applied after correcting our 
measured F606H' intensities for extinction and de-reddening the 
(C - Ti) co lours using F(C - Ti) = 1.97F(S - V) = 0.217 mag 
dHarris & C anterna 1979). 

We also estimate metallicities for our clusters from their 
[C-Ti) colour indices, using the relation o f lHarris & Harris! ( l2002h 
after correcting the colours for reddening: 



[Fe/H] = -6.037 [l - 0.82(C - FOo + 0.162(C - Fi)g] 



(5) 



This relation has been calibrated for classically old GCs, and it 
could give spurious metallicities for significantly younger clusters, 
if any such objects are in our sample. 

Table|4]lists the (V-F606)o colours and [Fe/H] v alues we have 
estimated in this way for the full ACS sample from IPaper III We 
have also added the 18 clusters from the sample of Paper I, which 
we have not re-observed. The first column of the table is the cluster 



© 2007 RAS, MNRAS 000,[TH26l 



D. 

Table 4. 



E. McLaughlin et al. 

Colours and photometric metallicities for 149 GCs in NGC 5128 " 



Name Tx (C-Ti) (C-r])o (V-F606)o [Fc/H] Name Ti (C-Ti) (C-r])o {V - F606)a [Fe/H] 

(1) (2) (3) (4) (5) (6) (1) (2) (3) (4) (5) (6) 



AATl 11563 


20.049 


1.091 


0.874 


± 


0.030 


0.091 


+ 


0.050 


-2.46 


+ 


0.1 


C146 


19.940 


1.773 


1.556 


± 


0.030 


0.177 


± 


0.050 


-0.70 


± 


0.1 


AAT113992 


19.818 


1.955 


1.738 


± 


0.030 


0.200 


± 


0.050 


-0.39 


± 


0.1 


C147 


20.055 


1.523 


1.306 


± 


0.030 


0.146 


± 


0.050 


-1.24 


± 


0.1 


AATl 15339 


19.561 


1.546 


1.329 


± 


0.030 


0.148 


± 


0.050 


-1.18 


± 


0.1 


C148 


20.207 


1.042 


0.825 


± 


0.030 


0.085 


± 


0.050 


-2.62 


± 


0.1 


AATl 17287 


20.450 


1.374 


1.157 


± 


0.030 


0.127 


± 


0.050 


-1.62 


± 


0.1 


C149 


19.680 


1.343 


1.126 


± 


0.030 


0.123 


± 


0.050 


-1.70 


± 


0.1 


AATl 18198 


19.031 


2.110 


1.893 


± 


0.030 


0.220 


± 


0.050 


-0.17 


± 


0.1 


C150 


19.780 


— 


1.400 


± 


0.300 


0.157 




0.063 


-1.00 


± 


0.6 


AATl 19508 


19.857 


1.832 


1.615 


± 


0.030 


0.185 


± 


0.050 


-0.59 


± 


0.1 


C151 


19.952 


2.317 


2.100 


± 


0.030 


0.246 


± 


0.050 


+0.05 


± 


0.1 


AATl 20336 


19.668 


1.809 


1.592 


± 


0.030 


0.182 


+ 


0.050 


-0.63 


+ 


0.1 


C152 


17.820 


— 


1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


AAT120976 


19.991 


1.511 


1.294 


± 


0.030 


0.144 


± 


0.050 


-1.27 


± 


0.1 


C153 


18.230 


— 


1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


C002 


17.937 


1.546 


1.329 


± 


0.030 








-1.18 


± 


0.1 


C154 


19.580 


— 


1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


C003 


17.081 


1.940 


1.723 


± 


0.030 


0.198 


± 


0.050 


-0.41 


± 


0.1 


C155 


21.358 


1.627 


1.410 


± 


0.030 


0.159 


± 


0.050 


-1.00 


± 


0.1 


C004 


17.498 


1.451 


1.234 


± 


0.030 


0.137 


± 


0.050 


-1.42 


± 


0.1 


C156 


17.651 


2.106 


1.889 


± 


0.030 


0.219 


± 


0.050 


-0.18 


± 


0.1 


C006 


16.510 


1.858 


1.641 


± 


0.030 


0.188 


± 


0.050 


-0.55 


± 


0.1 


C157 


19.210 


— 


1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


C007 


16.644 


1.534 


1.317 


± 


0.030 


0.147 


± 


0.050 


-1.21 


± 


0.1 


C158 


20.061 


1.606 


1.389 


± 


0.030 


0.156 


± 


0.050 


-1.05 


± 


0.1 


coil 


17.197 


2.011 


1.794 


± 


0.030 








-0.30 


± 


0.1 


C159 


19.929 


1.986 


1.769 


± 


0.030 


0.204 


± 


0.050 


-0.34 


± 


0.1 


C012 


17.358 


1.984 


1.767 


± 


0.030 


0.204 


+ 


0.050 


-0.34 


± 


0.1 


C160 


19.990 


— 


1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


COM 


17.407 


1.655 


1.438 


± 


0.030 


0.162 


+ 


0.050 


-0.94 


± 


0.1 


C161 


19.261 


1.953 


1.736 


± 


0.030 


0.200 


± 


0.050 


-0.39 


± 


0.1 


C017 


17.186 


1.422 


1.205 


± 


0.030 








-1.49 


± 


0.1 


C162 


20.909 


1.046 


0.829 


± 


0.030 


0.085 


± 


0.050 


-2.60 


± 


0.1 


C018 


16.891 


1.603 


1.386 


± 


0.030 


0.156 


± 


0.050 


-1.05 


± 


0.1 


C163 


20.431 


1.547 


1.330 


± 


0.030 


0.149 


± 


0.050 


-1.18 


± 


0.1 


C019 


17.554 


1.661 


1.444 


± 


0.030 


0.163 


± 


0.050 


-0.93 


± 


0.1 


C164 


19.600 


— 


1.400 


+ 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


C021 


17.397 


1.576 


1.359 


± 


0.030 








-1.12 


± 


0.1 


C165 


18.173 


1.931 


1.714 


± 


0.030 


0.197 


± 


0.050 


-0.42 


± 


0.1 


C022 


17.696 


1.516 


1.299 


± 


0.030 








-1.26 


± 


0.1 


C166 


20.720 


— 


1.400 


+ 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


C023 


16.686 


1.904 


1.687 


± 


0.030 








-0.47 


± 


0.1 


C167 


20.410 


— 


1.400 


+ 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


C025 


17.965 


1.952 


1.735 


± 


0.030 


0.200 


± 


0.050 


-0.39 


± 


0.1 


C168 


19.710 


— 


1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


C029 


17.533 


1.924 


1.707 


± 


0.030 


0.196 


± 


0.050 


-0.44 


± 


0.1 


C169 


20.389 


1.212 


0.995 


+ 


0.030 


0.106 


± 


0.050 


-2.08 


± 


0.1 


C030 


16.681 


1.789 


1.572 


± 


0.030 


0.179 


± 


0.050 


-0.67 


± 


0.1 


C170 


22.161 


1.148 


0.931 


± 


0.030 


0.098 


± 


0.050 


-2.27 


± 


0.1 


C031 


17.710 


2.023 


1.806 


± 


0.030 








-0.29 


± 


0.1 


C171 


20.670 


1.577 


1.360 


± 


0.030 


0.152 


± 


0.050 


-1.11 


± 


0.1 


C032 


17.854 


2.006 


1.789 


± 


0.030 


0.206 


+ 


0.050 


-0.31 


± 


0.1 


C172 


20.907 


1.242 


1.025 


± 


0.030 


0.110 


± 


0.050 


-1.99 


± 


0.1 


C036 


17.944 


1.378 


1.161 


± 


0.030 


0.127 


+ 


0.050 


-1.61 


± 


0.1 


C173 


21.057 


2.092 


1.875 


± 


0.030 


0.217 


± 


0.050 


-0.19 


± 


0.1 


C037 


17.962 


1.691 


1.474 


± 


0.030 


0.167 


± 


0.050 


-0.86 


± 


0.1 


C174 


21.424 


1.812 


1.595 


± 


0.030 


0.182 


± 


0.050 


-0.63 


± 


0.1 


C040 


18.490 


1.630 


1.413 




0.030 








-0.99 


± 


0.1 


C175 


21.838 


1.668 


1.451 


± 


0.030 


0.164 


± 


0.050 


-0.91 


± 


0.1 


C041 


17.969 


1.980 


1.763 




0.030 








-0.35 


± 


0.1 


C176 


21.297 


1.053 


0.836 


± 


0.030 


0.086 


± 


0.050 


-2.58 


± 


0.1 


C043 


18.068 


1.521 


1.304 


± 


0.030 


0.145 


± 


0.050 


-1.24 


± 


0.1 


C177 


21.087 


1.867 


1.650 


± 


0.030 


0.189 


± 


0.050 


-0.53 


± 


0.1 


C044 


18.148 


1.441 


1.224 


± 


0.030 








-1.44 


± 


0.1 


C178 


21.534 


1.481 


1.264 


± 


0.030 


0.140 


± 


0.050 


-1.34 


± 


0.1 


ClOO 


19.030 


1.409 


1.192 


± 


0.030 








-1.53 


± 


0.1 


C179 


21.043 


1.088 


0.871 


± 


0.030 


0.091 


± 


0.050 


-2.47 


± 


0.1 


ClOl 


20.077 


1.762 


1.545 


± 


0.030 








-0.72 


± 


0.1 


F1GC14 


19.671 


1.421 


1.204 


± 


0.030 


0.133 


± 


0.050 


-1.49 


± 


0.1 


C102 


21.038 


1.673 


1.456 


± 


0.030 








-0.90 


± 


0.1 


F1GC15 


19.531 


2.224 


2.007 


± 


0.030 


0.234 


± 


0.050 


-0.04 


± 


0.1 


CI 03 


18.439 


1.993 


1.776 


+ 


0.030 








-0.33 


± 


0.1 


F1GC20 


21.216 


1.839 


1.622 


± 


0.030 


0.185 


+ 


0.050 


-0.58 


± 


0.1 


C104 


19.403 


1.405 


1.188 


± 


0.030 


0.131 


+ 


0.050 


-1.54 


± 


0.1 


F1GC21 


21.358 


2.071 


1.854 


± 


0.030 


0.215 


± 


0.050 


-0.22 


± 


0.1 


C105 


21.758 


1.922 


1.705 


± 


0.030 


0.196 


± 


0.050 


-0.44 


± 


0.1 


F1GC34 


20.961 


1.896 


1.679 


± 


0.030 


0.193 


± 


0.050 


-0.48 


± 


0.1 


C106 


20.657 


1.879 


1.662 


± 


0.030 








-0.51 


± 


0.1 


F2GC14 


20.824 


1.678 


1.461 


± 


0.030 


0.165 


± 


0.050 


-0.89 


± 


0.1 


cm 


21.363 


1.172 


0.955 


± 


0.030 


0.101 


+ 


0.050 


-2.20 


+ 


0.1 


F2GC18 


21.147 


— 


1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


C112 


21.310 


1.668 


1.451 


± 


0.030 


0.164 


+ 


0.050 


-0.91 


+ 


0.1 


F2GC20 


21.013 


1.870 


1.653 


± 


0.030 


0.189 


± 


0.050 


-0.53 


± 


0.1 


CI 13 


19.120 


1.378 


1.161 


± 


0.030 


0.127 


± 


0.050 


-1.61 


± 


0.1 


F2GC28 


21.156 


1.701 


1.484 


± 


0.030 


0.168 


± 


0.050 


-0.84 


± 


0.1 


C114 


21.423 


2.037 


1.820 


± 


0.030 


0.210 


± 


0.050 


-0.27 


± 


0.1 


F2GC31 


20.154 


1.395 


1.178 


± 


0.030 


0.129 


± 


0.050 


-1.56 


± 


0.1 


C115 


19.507 


1.372 


1.155 


± 


0.030 


0.127 


± 


0.050 


-1.62 


± 


0.1 


F2GC69 


19.301 


1.813 


1.596 


± 


0.030 


0.182 


± 


0.050 


-0.63 


± 


0.1 


C116 


21.837 


1.967 


1.750 


± 


0.030 


0.202 


± 


0.050 


-0.37 


± 


0.1 


F2GC70 


20.039 


1.253 


1.036 


± 


0.030 


0.112 


± 


0.050 


-1.96 


± 


0.1 


C117 


19.262 


2.011 


1.794 


± 


0.030 


0.207 


± 


0.050 


-0.30 


± 


0.1 


G019 


18.636 


1.422 


1.205 


± 


0.030 








-1.49 


± 


0.1 


C118 


20.672 


1.924 


1.707 


± 


0.030 


0.196 


± 


0.050 


-0.44 


± 


0.1 


G170 


18.727 


1.844 


1.627 


± 


0.030 


0.186 


± 


0.050 


-0.57 


± 


0.1 


CI 19 


20.669 


0.487 


0.270 


± 


0.030 


0.015 


± 


0.050 


-4.77 


± 


0.1 


G221 


18.828 


1.715 


1.498 


± 


0.030 


0.170 


± 


0.050 


-0.82 


± 


0.1 


C120 


21.433 


1.828 


1.611 


± 


0.030 


0.184 


± 


0.050 


-0.60 


± 


0.1 


G277 


18.607 


1.530 


1.313 


± 


0.030 








-1.22 


± 


0.1 


C121 


22.446 


1.080 


0.863 


± 


0.030 


0.090 


± 


0.050 


-2.49 


± 


0.1 


G284 


19.412 


1.853 


1.636 


± 


0.030 


0.187 


+ 


0.050 


-0.56 


± 


0.1 


C122 


22.330 


— 


1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


G293 


18.693 


1.349 


1.132 


± 


0.030 


0.124 


+ 


0.050 


-1.69 


± 


0.1 


C123 


20.244 


1.637 


1.420 


± 


0.030 


0.160 


± 


0.050 


-0.98 


± 


0.1 


G302 


18.728 


1.450 


1.233 


± 


0.030 








-1.42 


± 


0.1 


C124 


21.574 


1.498 


1.281 


± 


0.030 


0.142 


± 


0.050 


-1.30 


± 


0.1 


K131 


18.741 


2.104 


1.887 


± 


0.030 


0.219 


± 


0.050 


-0.18 


± 


0.1 


C125 


20.664 


1.671 


1.454 


± 


0.030 


0.164 


± 


0.050 


-0.91 


± 


0.1 


PFFOU 


18.553 


1.398 


1.181 


± 


0.030 


0.130 


± 


0.050 


-1.55 


± 


0.1 


C126 


22.663 


1.354 


1.137 


± 


0.030 


0.124 


± 


0.050 


-1.67 


± 


0.1 


PFF016 


19.334 


1.853 


1.636 


± 


0.030 


0.187 


± 


0.050 


-0.56 


± 


0.1 


C127 


21.620 


1.593 


1.376 


± 


0.030 


0.154 


± 


0.050 


-1.08 


+ 


0.1 


PFF021 


18.814 


1.318 


1.101 


± 


0.030 


0.120 


± 


0.050 


-1.77 


± 


0.1 


C128 


20.947 


2.131 


1.914 


± 


0.030 


0.222 


± 


0.050 


-0.14 


± 


0.1 


PFF023 


18.971 


1.511 


1.294 


± 


0.030 


0.144 


± 


0.050 


-1.27 


± 


0.1 


C129 


20.831 


1.757 


1.540 


± 


0.030 


0.175 


± 


0.050 


-0.73 


± 


0.1 


PFF029 


19.098 


0.621 


0.404 


± 


0.030 


0.032 




0.050 


-4.20 


± 


0.1 


C130 


19.851 


1.368 


1.151 




0.030 


0.126 


± 


0.050 


-1.63 


± 


0.1 


PFF031 


18.952 


1.479 


1.262 


± 


0.030 


0.140 


± 


0.050 


-1.35 


± 


0.1 


C131 


20.202 


1.363 


1.146 


± 


0.030 


0.125 


± 


0.050 


-1.65 


± 


0.1 


PFF034 


19.384 


1.502 


1.285 


± 


0.030 


0.143 


+ 


0.050 


-1.29 


± 


0.1 


C132 


18.896 


1.480 


1.263 


± 


0.030 


0.140 


+ 


0.050 


-1.34 


± 


0.1 


PFF035 


19.132 


2.000 


1.783 


± 


0.030 


0.206 


± 


0.050 


-0.32 


± 


0.1 


C133 


19.300 




1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


PFF041 


19.007 


1.471 


1.254 


± 


0.030 


0.139 


± 


0.050 


-1.37 


± 


0.1 


CI 34 


20.711 


1.735 


1.518 


+ 


0.030 


0.172 


± 


0.050 


-0.78 


± 


0.1 


PFF052 


19.416 


1.442 


1.225 


± 


0.030 


0.135 


± 


0.050 


-1.44 


± 


0.1 


C135 


18.900 




1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


PFF059 


19.479 


1.943 


1.726 


± 


0.030 


0.199 


± 


0.050 


-0.41 


± 


0.1 


C136 


21.190 


1.295 


1.078 


± 


0.030 


0.117 


+ 


0.050 


-1.84 


+ 


0.1 


PFF063 


19.452 


1.176 


0.959 


± 


0.030 


0.102 


± 


0.050 


-2.19 


± 


0.1 


C137 


19.040 




1.400 


± 


0.300 


0.157 


+ 


0.063 


-1.00 


+ 


0.6 


PFF066 


19.436 


1.629 


1.412 


± 


0.030 


0.159 


± 


0.050 


-1.00 


± 


0.1 


C138 


19.974 


1.555 


1.338 


± 


0.030 


0.150 


± 


0.050 


-1.16 


± 


0.1 


PFF079 


19.140 


1.484 


1.267 


± 


0.030 


0.141 


± 


0.050 


-1.33 


± 


0.1 


C139 


18.863 


1.866 


1.649 


± 


0.030 


0.189 


± 


0.050 


-0.53 


± 


0.1 


PFF083 


19.436 


1.683 


1.466 


± 


0.030 


0.166 


± 


0.050 


-0.88 


± 


0.1 


C140 


19.830 


1.570 


1.353 


± 


0.030 


0.152 


± 


0.050 


-1.13 


± 


0.1 


R203 


19.384 


1.240 


1.023 


± 


0.030 


0.110 


± 


0.050 


-2.00 


± 


0.1 


C141 


20.945 


1.288 


1.071 


± 


0.030 


0.116 


+ 


0.050 


-1.86 


± 


0.1 


R223 


18.186 


1.708 


1.491 


± 


0.030 


0.169 


± 


0.050 


-0.83 


± 


0.1 


C142 


17.640 




1.400 


± 


0.300 


0.157 


± 


0.063 


-1.00 


± 


0.6 


WHH09 


18.300 


1.979 


1.762 


± 


0.030 


0.203 


± 


0.050 


-0.35 


± 


0.1 


C143 


20.525 


1.345 


1.128 


± 


0.030 


0.123 


± 


0.050 


-1.70 


± 


0.1 


WHH16 


18.599 


2.036 


1.819 


± 


0.030 


0.210 


± 


0.050 


-0.27 


± 


0.1 


C144 


21.959 


1.257 


1.040 


± 


0.030 


0.112 


± 


0.050 


-1.95 


± 


0.1 


WHH22 


18.032 


1.631 


1.414 


± 


0.030 


0.159 


± 


0.050 


-0.99 


± 


0.1 


C145 


17.814 


1.510 


1.293 




0.030 


0.144 


± 


0.050 


-1.27 


± 


0.1 



























" A blank entry in Column (5) indicates a cluster appearing in the sample of |Harris et alj |20Q2| | but not observed by us here, and for which we therefore do not need a { V - f 606)o colour. 
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Figure 2. V-band mass-to-light ratio vs. metallicity, from the population- 
synthesis models of Bruzual & Chai'lol (2003) for single-burst clusters with 
the stellar IMF of Chabrier 12003). Curves are for cluster ages of 1,9, 11, 
13, and 15 Gyr. Note that all of our clusters have [Fe/H] < 0, and thus 
T^"'' < 4 Mq Lq'j/ for an age of 13 Gyr. A typical value is closer to T^"*" ^ 
2.1 Mo for an average <[Fe/H]> ^ -1.0 in NGC 5128. 



ID. The second and third columns are taken directly from Tables 1 
and 2 of Paper II; they are the observed Washington T\ magnitudes 
and {C -Ti) colours (not corrected for extinction/reddening). The 
fourth column is the de-reddened (C - ri)o = (C - Ti) - 0.217. We 
have adopted a uniform uncertainty of ±0.03 mag on this colour 
for most clusters. Note that there are 16 clusters (all in our sub- 
set of newly discovered objects) for which there is no observed 
{C —Ti) index. To these 16 objects we have assigned (C - ri)o = 
1.4 ± 0.3 mag, which is the average and dispersion of the measured 
(C - rOo for the other clusters. Columns 5 then gives the colour 
(V - F606)o; it is left blank for the 18 clusters from Paper I, whose 
surface-brightness profiles we are not modeling here. Column (6) 
reports the [Fe/H] inferred from equations Q and ([5j- 

There are two objects (CI 19, PFF029) with colours so blue 
that equation Jsj implies [Fe/H] < -4. These are likely to be much 
younger objects than the classical ~ I3-Gyr globular cluster ages 
(see Paper II), and their metallicities should not be taken seriously. 



3 MODELS 

3.1 Population synthesis models: MIL ratios 

Velocity dispersions have been measured for 20 of the GCs in 
the combined sample of this paper and IPaper ll l lMartini & Hoi 
l2004lReikuba et alj|2007h . These can be used in conjunction with 
surface-brightness modeling to estimate dynamical mass-to-light 
ratios. However, since this is possible for only for a small minority 
of our clusters, we have chosen instead to use population-synthesis 
models to predict mass-to-light ratios for the full sample, and then 
to use these to produce (for example) predicted velocity disper- 
sions that can be compared to current and future spectroscopic data. 
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Figure 3. Histogram of population-synthesis V-band mass-to-light ratios 
for 145 GCs in NGC 5128 (excluding C145, C152, C156, and C177; see 
Table |2j, for an assumed age of 1 3 Gyr; and the same distribution for 
148 Galactic globulars (normalized to a total population of 145), from 
iMcLaughlin & van der Marell f20o3) . 



Since we have already derived the requisite corrections to put our 
photometry on the standard V scale, we need only calculate the 

mass-to-light ratios in this bandpass. 

This approach is also taken bylM cLaughlin & van der Marell 

( l2005h in their study of 153 globular and young massive clusters 
in the Milky Way, the Larg e and Small Magellanic Clou ds, and 
the Fornax dwarf spheroidal. iMcLaughlin & van der Marell discuss 
the differences in population-syn thesis mass-to-light ratios , Ty"*", 
obtained using different models ( iBruzual & Charlol |2003| versus 
iFioc & Rocca- Volmerange|[l997l) and for a number of different as- 
sumed stel lar IMFs. They ultimately work with the results from 
the code of iBnizual & Chariot! ( I2OO3I) using the (disk-star) IMF of 
IChabriej l l2003r to produce catalogues of dynamical properties for 
all of their clusters. We thus do the same here. 

Figure|2]illustrates model curves of T^"'' versus cluster metal- 
licity, for ages from 7 to 15 Gyr. Given an age, and a metallicity 
from Table (4] we interpolate on the curves in Figure |2| to obtain a 
model T^"'' for any cluster in our sample. We eventually assume an 
age of 13 Gyr for all GCs, in ^and subsequent analyses. 

Table|5]repor ts Yj'"'' v alues for our ACS sample plus the 18 ad- 
ditional GCs from jPaper ll The first column lists the cluster name, 
and subsequent columns give the model mass-to-light ratios for as- 
sumed ages of 7, 9, II, 13, and 15 Gyr. As was just suggested, in 
the majority of cases we have assumed the cluster metallicity in 
Column (6) of Table [4] to evaluate T'^"''. However, there are 9 clus- 
ters in the tab le with apparent [Fe /H] < -2.25, below the minimum 
for which the lBruzual & CharloJ models are defined. We have reset 
these to have [Fe/H] = -2.25 so we can still obtain some estimate 
of T'^"''; but note that such low metallicities follow from very blue 
(C - Ti)o colours, which may be indicating that these objects are 
much younger than 13 Gyr (CI 19 and PFF029, which we men- 
tioned at the end of 92.31 are in this group). 

Figure [3] shows a histogram of t^y'^ for GCs in NGC 5128 as- 
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Table 5. Model V-band mass-to-light ratios for 149 GCs in NGC 5128 



Name 



7Gyr 



9Gyr 



11 Gyr 



13 Gyr 



15 Gyr 



Name 



7 Gyr 



9 Gyr 



11 Gyr 



13 Gyr 



15 Gyr 



AATl 11563 " 
AATl 13992 
AATl 15339 
AATl 17287 
AATl 18198 
AATl 19508 
AAT120336 
AAT120976 
C002 
C003 
C004 
C006 
C007 

coil 

C012 

COM 

C017 

C018 

C019 

C02I 

C022 

C023 

C025 

C029 

C030 

C031 

C032 

C036 

C037 

C040 

C04i 

C043 

C044 

ClOO 

ClOl 

C102 

C103 

C104 

C105 

C106 

Clll 

C112 

C113 

C114 

C115 

C116 

C117 

C118 

CI 19° 

C120 

C121 ' 

C122 

C123 

C124 

C125 

C126 

C127 

C128 

C129 

C130 

C131 

C132 

C133 

C134 

CBS 

CI 36 

CI37 

C138 

C139 

C140 

C141 

C142 

C143 

C144 

CI45 



1 18"*"°'^^ 

I 7g+8:?i 

'■'^-0.12 
1 9^+0.03 
-0.03 

1.57+0-1/, 



1.53" 



-0!09 



1 91+0.03 
^■^ -0.02 



1 9-1+0.03 



-0!03 
, 77+0.12 
^■"-0.12 
I in+0.02 
^■^-^-0.01 

1.61+0-'2 
1 89+^-^- 

° -0.12 

^ 85+0.12 

, T,+8:i^ 

-0.05 
1 iq+0.01 

27+0-05 

1 09+8:8^ 

^■^^-0.05 

1.25+004 



1.21 + 
-u.uz 

i.7o;»!f 



1 79 . 

-0.12 
1 74+0.12 

1 so+'^-l^ 

-0.09 
1.91+"'5 

l-8»-0.12 
I lg+0.01 

J 3j+8:8'> 
1 9g+8:8§ 

'■^^-0.04 
J g4+8.12 

1 -,i+8:A? 

-0.02 
[ jj+0.02 

J 18+8:81 
1 A'i+SSi 

-^-0.08 
1 33+0.07 

1.86 ,, ,2 
I 1 0+0.01 

'■'=-0.01 
I 74+''12 

S6;8:13 

l.OD 

1 18+0.01 
1 ^^+8:81 

'■^^-0.05 
1 18"^"'" 
1 q3+8:!J 
'•"-0.12 
1 18+''-'" 
1 01+8:?^ 
1-81-0.12 

1 89+'"2 
'-°^-0.12 
1 74+''-'2 
'-' -0.12 
1 lfi+0.01 
'-'"-O.Ol 

'■-'°-o.io 

1 1 0+0.01 
'-'=-0.01 
1 o()+0.50 
'-^'-0.11 

'-^"-0.04 
12o+(l-»2 

1 T^+8:8^ 

'-^^-0.05 

J Ig-*-"-"' 

1 26+'^-'^^ 
'-^°-0.03 
207+0.13 

] 45+8:45 

1 10+8:8? 

1-18-0.01 

1 1 0+0.01 

'-'=-0.01 
, 20+0-02 

1 o()+8:SS 
'-^'-0.11 

1 41+0.08 
'-^'-0.07 
1 2Q+0.50 
'-^'-0.11 
1 10+0.01 
'-'8-0.01 
1 29+0.50 

1 ')^+8:il 

'-^^-0.03 
1.64+0-12 

, 24+8:il 

'- -0.03 

1 10+0.01 

'-'8-0.01 
1 29+0.50 
'-''^-O.ll 
1 10+0.01 
'-'8-0.01 
1 10+0.01 
1-18-0.01 
1 21+0.03 
'-^'-0.02 



145+0.01 

T i7+8:?i 

^-"-013 
1 50+0-05 
'--^ -0.04 

1 42+0.01 

2 4c+8:?i 

^-^^-0.13 
1 92+013 

, oo+8:li 

1-88-0.11 

1.47+r"' 



7+0:04 

'-0.03 

1 50+0-05 
'--^ -0.04 

2 14+013 
•^-'+-0.13 
I 44+0.02 
'- -0.02 

1 97+013 

1 49+8:41 

T '10+8:?^ 
2.28 „ 

2 23+0- 13 

1 61+''-''^ 
'-"'-0.06 
I 43+0.02 
'-*^-0.01 
1 gg+0.06 



:8:15 



1 £:7+0.08 
1 52+0.05 

1 48+8:8^ 
'■+°-0.03 
207+0.13 

2 17+8:1^ 
^■"-0.13 
2 11+0.13 

1 84+'^-! 5 
'■ -0.10 

2.30 

2.27_ 

1 42+0:61 

1 57+8:84 

, cQ+8:8? 

'-■'8-0.06 

2 22+0.13 

1 

'■+"-0.03 
1 44+0.02 

1 43+8:85 
1 7o+8:?l 

'■'8-0.09 

1 64+0-OS 
'-°*-0.07 

2 2^+0.13 
^-"-0.13 
I 43+0.02 
'-+^-0.01 
2.11 



+0.13 

Ml 



1 42+0.01 

2 20+^-'i 
^-^"-0.13 
2 29+0- 13 
^-^'-0.13 
2 ll+O '' 
^-"-0.13 
I 46+0 01 
'■™-0.02 
191+0.13 
'^'-0.11 
1 46+0^0' 
'■^"-0.02 

1.58 
1.59+ 
1.46+ 
1.63+ 



+0.58 
-0.16 



1.78_ 

1 42+S8I 
'■+^-0.01 
1 42+0-01 
'-^^-0.01 
145+0.03 

1 ^o+8:§i 

1-58-0.16 
I 7^+0.10 
'-'^-0.08 
1 ^0+0.58 
1-58-0.16 
1 43+0.01 
'-^^-0.01 



1.58 

i.5r 



+0.58 



0.04 
Q+0.13 



1.52+ 



l-''9*.., 

:45 

^-0.04 
I 43+0.01 
'-*^-0.01 
1 ^8+0.58 
1-58-0.16 
I 42+0.01 
'-^^-0.01 
I 43+0.01 
'-^^-0.01 
I 47+0.04 
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1 59+0.01 

T r^+8:?^ 

■^-^^-015 
1 77+0.05 
'-"-0.04 

1 57+0.02 

2 89+''-''^ 
7 25+014 

2 20+8:1^ 
-012 
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'-'^-0 03 



i5 



1 76+0-"- 
'-' -0.04 
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1 69+'^-'''^ 
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230+0.15 

1 75+8:4^ 
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1 59+0.03 
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'■"'-001 

2.10+1 



1.93;! 



-Q.IO 
, ,-0.09 
- -^-0.08 

2 63+0" 
^■"^-0.16 

1 67+002 
'■°'-0.01 
245+0.16 

2 36+^-'^ 
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' -0.12 
1 0+0.01 
^*^-0.01 
23+0.03 

9Q+8:§^ 

^^-0.11 
.76+0-12 

29+8:^^ 



-0.12 
,+0.50 



?8:1? 



54+0.12 
"^-0.11 
20+0.02 

18+8:8? 

19+8:81 
-0.01 

19+014 
'-^-013 
58+011 
JO n in 



a;8:li 



53+0.11 
1 8+8:8? 

'8-0.01 
19+0.01 
"-0.01 
59+0" 



-0.02 
0+0.11 
-0.10 



18+0.01 
'8-0.01 
19+0.02 
'-"-0.01 
2 03+0.12 
■ -Q.12 
1.18+^ 



"-0.01 



n+0.1 



21+0.03 

18+8:8? 

70+8:84 

^"-0.02 
7 1 +0.03 

^'-0.02 
87+0.12 

70+8:43 

^"-0.01 
19+0.02 
"-0.01 
77+0.12 
' '-0.12 
18+0.01 
'8-0.01 
29+0.05 

7n+8:8l 
^" -0.02 



+0:07 



38+0.08 

84+8:'^ 
"'-0.12 
93+0.12 
"-0.12 
79+0.05 
-0.04 



1 81+0" 

1 48+8:^3 
'■*°-0.03 
1 46+00' 
'■+"-0.02 

1 42+0.01 

1 58+8:5^ 

2 80+8:^^ 
I ^8+8:^^ 

'■58-0.16 
1 58+0.58 

I .^8+8:^^ 

158-0.16 

1 58+0.07 

2 44+8:5? 
^■++-0.13 

1 58+0.58 

1 55+8:4^ 

2 77+8:5? 

^■"-0.13 

1 58+0.58 

2 17+8:15 
^■' '-0.13 

1 46+00' 
'■^°-0.02 
150+0.05 

, rD+8:§i 

1-58-0.16 

2 13+0.13 

1 ^o+8:4i 

1-58-0.16 
1 58+0-58 

1 ^o+8:i^ 

1-58-0.16 
1 44+0.01 

1 45?8:81 

'-+"-0.02 
153+0.05 

1 44+8:8! 

'-++-0.01 
243+0.14 

1 88+8:15 
1 ^7+8:44 

1.63 ,j„7 
145+0.01 

1 99+8:53 

'-"-0.12 

1 45+0.03 

145+8:8? 
1 43+8:83 

'-^^-0.01 

2 54+0.18 

1 93+8:15 
'- -an 

7 39+0.14 
-0.13 
2.06+O'3 

1 65+8:''^ 

'-"-'-O.OT 
1 58+0.58 

1 99+8:15 

'- -012 
1 58+0.09 
1.68 

1 43+0.81 
'-'^-0.01 
1 88+0-12 

1 43+8:41 
'-+^-0.01 
1 43+0.02 
'-*^-0.01 



1 95+0:13 
'-"-0.12 
1 70+0-0!> 
'-'"-0.08 
1 49+0.04 
■ -0.03 
1 96+0' 3 
'■'"-0.12 



:83 



1 42+0.01 
'■+^-0.01 
1 44+0.02 
'■++-0.02 
244+0.15 

1 43+8:4? 
'■+^-0.01 
1.96;0'3 

1 42+8:4? 
'■+^-0.01 
1 47+0.04 

■ -0.03 
1.46+ 

1 45+-.- 

■ -0.02 
1 47+0.03 
'■^'-0.03 
225+0.13 

1 45+8:4? 

■ -0.02 

1 44+0.02 

■ -0.02 

2 14+0'3 
^■'+-0.13 
1 45+0.01 
'■^■^-0.02 
1 58+0.07 

1 46+8:85 

'■+"-0.02 
,55+0.08 

I 44+8:8? 
'■++-0.01 

1 69+009 
' ■°'-0.08 

2 72+0.13 

-0.13 
2 72+0.13 
^■^^-0.13 
I C8+0.07 
158-0.06 



2 12+0.12 

1 74+8:44 

'■'*-0.04 
1 69+00' 
'■"'-0.02 
1.66+?^9| 



'-0.20 



2.87 „;j7 

1 83+0.07 

2 61+8:?^ 

-0.16 
J gg+0.66 

9 cT+8!?8 
^■■^^-0.15 
1 69+0-01 
'■'^^-0.02 
I y-^+0.05 

1 Q(;+8!^^ 

l.86_o_2o 



+0.16 



2.49 
1 
1 
1 

1 58+001 

1 69+8:8! 

'■"'-0.02 
I 80+0-06 
1-80 „.g5 

' -0.01 



1.67+' 



+0.19 



2.85_ 
2.20+j 
1.93+! 



9 To+0.15 
^■^^-0.13 
J y J +0.04 

^ 69+8:8? 

^■^^-0.02 
3 14+0-23 
-0.20 

2.26:0» 

2.80+0;|? 
240+O15 

1 94+8:45 
-0.08 
1 85+0.66 

7 77+8:?5 

^-^ -013 



-0.01 

1 58+0.02 
'-°8-0.02 
7 78+0.14 

2.28 „ ,3 

2 OO+O-'O 

-0.09 

1 7^+0.05 
'■'■^-0.04 

2 79+0.14 
^■^'-0.13 
1 66+00' 
'■™-0.01 
1 69+003 
'■°'-0.02 



+0:19 



2.8T 



'■'5«:!!:ol 

1 73+0.04 

1 69+8:8^ 

1.69 „g, 

'-0.03 
1 77+0.04 
'■'^-0.03 



1.71 + 



2.64;0'7 

1 7o+8:i'i5 

■ -0.03 

1 69+003 
'■"'-0.02 

2 50+0 ' 6 
^■■'"-0.15 
1 69+00' 
'■"'-0.02 



+0.08 



1.86 



■ ' -0.01 
I 99+0.10 
'■"-0.09 



244+0.14 

I 99+8:4i 

'■"-0.05 
I 94+0.01 
'■'^-0.02 



1 

2 i3+i):W 

3 91+8:35 
2 13+8:^^ 

^■'^-0.25 
2 13+0.77 

2 13+8:3? 
2 1 7+8:55 

■ -0.08 
7 77+0.24 
^■^^-0.22 
2 13+0.77 

2 10+8:1 

3 01+8:5^ 

^■"'-0.19 
2 13+0.77 

2 91+8:35 

^■"-0.18 
1 94+0.01 
'■'+-0.02 
202+0-07 

7 17+8:5? 

^-'^-0.25 
2.86+9-1? 



■?8:« 

-0.25 
2.13+9-" 



2.13+ 



2.13+1 



-0.25 
1 91+0.02 

I 94+8:8? 
'-'+-0.02 
205+0-08 

I 90+8:85 

'-'"-0.02 
3 3 1 +0.24 



?8:?4 



2.53_ 
7 7i+y:t5 

^-^'-0.09 
1 94+0.01 

7 *;o+8:53 
2.68ui5 

1 94+0.05 



1.94; 
1.90 



:8:83 
-0.02 



3.67+028 

2 60+8:16 
^■""-0.14 
7 7C+O.23 
^■^^-0.21 
275+018 

7 23+8:15 
-0.09 
2 13+0.77 

7 58+8:?? 
2.68 „ 15 

2 28+0.12 

1 89+8:i'i5 

'-"'-0.01 
2.53+0-15 



1 90+' — 
'-'"-0.01 
1.90+0-03 



'"-655 



2 62+o:'6 

-015 

2 30+012 
^-^"-0.10 
1 99+0.06 
'■"-0.05 
253+0.17 

1 88+8:4? 

'■88-0.01 
1 92+0.04 
'■'^-0.03 

3 33+0.24 

1 89+8:23 

'■"'-0.02 
253+0.17 

1 88+8:4? 

188-0.01 
1 97+0.05 

1 94+8:8? 
1 94+8:81 

■ -0.03 
1 96+005 
'■'"-0.04 
305+0.21 

1 93+8:45 

■ -0.03 
1 92+0.04 

■ -0.03 
7 88+0.19 
^■88-0.17 

1 93+0.02 
'■'^-0.02 

2 13+0.09 

1 95;8:8§ 

'■'■^-0.04 
2.24+0" 

1 90+8:8^ 

'■'"-0.02 
2 79+0.12 
^■^'-0.10 

2 99+0.20 
^■"-0.18 

3 14+0.22 
^■'^-0.20 
2 14+0.09 

■ -0.08 



2 72+0-17 

7 7 1 +8:4^ 

^-^'-0.04 
7 1 8+0.01 
'^-' 8-0.02 
2 13+0.01 

7 77+8:84 

^-^'-0.24 
437+0.34 

2 37+8:55 

7 77+8:5i 

^-^'-0.24 
2 77+0.92 
^-^'-0.24 
2 77+0.10 
^-■"-0.08 
3.76+0| 

2 33+0.09 

3 40+8:5^ 
^-+"-0.22 

2 37+0.92 

3 30+8:^^ 

^-^"-0.21 
7 18+0.01 
^-'8-0.02 
224+0-07 

7 77+8:8? 

^-^'-0-24 
3 23+0-22 

7 77+8:54 

'■■^ '-P-24 
7 77+0-92 
■^■^'-0-24 
2-37* 



-0-24 
+0-02 



2-16" 



2-83+ 
2-45+ 



2 18+0.04 

2 18+8:8? 
2 14+8:83 

^■'+-0-02 
4 11+0-30 

7 97+8:38 

^■' -0-18 

3 67+0 24 
^■"'-0-23 



i+B:5T 

-0-19 



2-14+1 
2-94+] 
2-55+1 

2-22+' 



.12 
.06 
^-0.04 
2 95+0.20 

^■'■'-o.is 

2 17+0.01 
^■'^-0.01 
2.16;r" 
3.76 



ft+o:o3 

0.02 
0.25 

:8:S? 



220+U-05 

2 18+8:8? 

7 17+8:84 
^-"-0.03 
2 19+0.05 



m 



2-17_ 

7 ir+(^:(i3 

^-'-^-0.02 
3 26+0-22 
^-^"-0.21 
7 17+0.01 
^-' '-0.02 
2 37+0.10 

7 1 8+8:81 

248+0.13 

2 1 c+8:45 

^■'■'-o.oi 

2 54+0.14 
^■^*-0.11 
7 78+0.22 
j.3S 

3 56+0 23 
^■^"-0.22 
2 37+0.10 

■ ' -0.08 



" Denotes clusters inferred on the basis of (C - T\ )o colour to have [Fe/H] < -2.25 in TableRl Mass-to-hght ratios are computed for them assuming [Fe/H] = -2.25 instead- 



© 2007 RAS, MNRAS 000,[TJ{26] 



Globular clusters in NGC 5128 9 



suming a common age of 13 Gyr, and compares it to the distribution 
of mass-to-light ratios for 148 Galac tic globulars, as calculated by 
iMcLaugh lin & van der M aref (2005). Note that the distribution in 
NGC 5 128 peaks at a slightly higher T!^"'' and is rather broader than 
that in the Milky Way. This just reflects the dilferent averages and 
dispersions of the GC metallicity distributions in the two galaxies. 

3.2 Structural models 

We fit a number of models to the density profile of each cluster 

observed with ACS/WFC 

First is the usua l lKingHl966l) single-mass, isotropic, modified 
isothermal sphere, which is defined by the stellar distribution func- 
tion (phase-space density) 



( expi-E/a-g] - 1 , £ < 
\0 , £>0, 



(6) 



where E is the stellar energy. Under certain conditions, this formula 
roughly approximates a steady-state solution of the Fokker-Planck 
equation (e.g.. King 1965). It is the standard model that is routinely 
fit to GC surface-brightness profiles. 

Second is a further modification of a single-mass, isotropic 
isothermal sp here, based on a model originally introduced by 
IWilsonI ( Il975h for elliptical galaxies. It has recently been fitted to 
globula r and young massive clusters in the M ilky Way and its satel- 
lit es bvlMcLaughlin & van der Marell ilOO^ . and to GCs in M31 
bv lBarmbv et alj ( |2007|) ." This model is again defined by a phase- 
space distribution function: 



( exp[-£/o-2] - 1 + E/crl , £ < 



(7) 



\0 , £ > 

IWilsonI ( [1973) originally included a multiplicative term in the dis- 
tribution function depending on the angular momentum J-, in order 
to create axisymmetric model galaxies. We omit this term to make 
spherical and isotropic cluster models, but we st ill refer to e quation 
(|7J as Wilson's model. The connection with the lKing| ( ll966h model 
in equation (|6j is clear: the extra +E/(Tg in the first line of equation 
(|7j removes the linear term from a Taylor series expansion of the 
Maxwellian exp{-E/(T^^) near the zero-energy (tidal) boundary of 
the cluster. Its effect is to make clus ters that are spatially more ex- 
tended than in isotropic iKind d 19661) models. To see this, note that 
adding E/a^ to the exponential term causes a much more signifi- 
cant decrease in the phase-space density of tightly bound stars with 
£ <K than it does in the density of loosely bound stars with E near 

0, whic h have large orbital apocentres. 

A s [ McLaughlin & van der Marell ( l2005h discuss in detail, ex- 
tended [wilsonl -tvpe cluster models fit the majority of old and young 
massive clusters in the Local Group at least as well as, and of- 
ten significantly better than. King models. The difference between 
the King and Wilson models when compared with real cluster 
data tends to show up most prominently at large radii approaching 
the nominal tidal radius. At small and intermediate radii (within 
the core, and out to slightly beyond the half-mass radius), the 
two model profiles are quite simila r. It is worth noting (see also 
[McLaughlin & van der Marelll2005l) that in the earlier era of mea- 
sured data that almost never extended to large radii, it would not 
have been possible to decide between these models. The much 
larger amounts of cluster profile data now available, which ex- 
tend more reliably out to the cluster tidal radii and even beyond, 
make it possible to start discriminating betwen models more accu- 
rately. In such cases, we find that the Wilson model treats the en- 
tire cluster profile more accurately without the need to invoke arbi- 



trary amounts of "extra-tidal light" beyond the formal King profile 
boundaries. 

Spherically symmetric, dimensionless densit y and v elocity- 
dispersion profiles are obtained for King ( 1 96^ and ! Wilson[ models 
by appropriate integrals of f(E) over all velocities in the first case, 
and over all spatial radii in the second. These are then integrated 
along the line of sight to produce normalized surface densities I/Iq 
and velocity dispersions CTp/cro as functions of a dimensionless pro- 
jected clustercentric radius R/rg, for comparison with observations. 
Here fq is a scalelength associated with, but not equivalent to, the 
observed half-power radius R^, as discussed below. The shapes of 
both profiles are fully specified by the value of a dimensionless 
central potential. Wo = -(f)(0)/o-Q > 0. In principle Wq can take 
on any real value between and oo, with the latter limit corre- 
sponding in both models to a non-singular isothermal sphere of 
infinite extent. Wo bears a one-to-one relationship with the more 
intuitive concentration parameter: c = log (f,/ro), where r, is the 
tidal radius of the model cluster [p(r,) = 0]. As suggested above, 
however, the tidal radii of iWilsoni mode ls are generally larger than 
those of otherwise similar Kind ( [l963) models, so that the same 
cluster data will almost always return different c values for the two 
models; this parameter, which is frequently mentioned in discus- 
sions of GC properties, is a h i ghly m odel-dependent quantity. See 



I, IS a tugl 

[m cLaughlin & van der Marell I2OO5I) for further discussion of this 
point. 

The third model we fit to our data is defined by directly 
param etrizing the observable surface -density profile. It is the Se rsid 
jl968[) or model, which, like the lWilsonI model, was also orig- 
inally designed for application to galaxies. We write this slightly 
differently than is usually done: 



I(R) = Io exp[-ln(2)x(i?/ro)'/"l 



(8) 



in which n > 0, and rg is the projected radius a t which the surface 
density falls to half its central value Iq. Often 'Sersic models are 
defined in terms of the projected half-light (or effective) radius — 
which we denote by Ri, — such that the exponentiated term in brack- 
ets in equation (|8j is [-b„(R/R,,)^'"], with b„ = ln(2) xj^/ro)"" a 
function of n that we compute numerically. Note that lSersi^ models 
have a formally infinite spatial extent. However, the density profile 
falls steeply at very large R, and the models thus have finite total 
luminosities and well-defined Ri, for any n. One slight complica- 
tion is that the de-projected density profiles j(r) corresponding to 
the observable I(R) are weakly divergent in the central limit r — > 
for « > 1. We return to this point just below. 

As we hav e already discussed, setting Wo at some value for 
lKind ( ll966l) and lWilso nl models leads to the full definition of three- 
dimensional and projected density and velocity-dispersion profiles, 
the shapes of which are fixed by Wq or, equivalently, a va lue of 
c = log(r,/ro). The analogous shape parameter for the [Sersi^ mod- 
els is the index n in equation (|8j. Even though there is no longer 
a connection between the shape parameter and any tidal radius in 
this case, there is still a one-to-one correspondence between n and 
a dimensionless central potential Wo = -0(O)/o'o, which is always 
finite. This works in the sense that a higher n implies density pro- 
files that fall off more slowly with increasing radius, which in turn 
implies a deeper central potential, or higher Wo. 

Given a value for n, we numerically deproject equation ^ 
to compute the volume luminosity density /(r), which is then 
used to solve the spherical leans equation ( jsinnev & Tremain^ 
Il987h — assuming unit mass-to-light ratio and an isotropic velocity 
ellipsoid — to obtain a normalized velocity-dispersion profile. This 
can be re-projected along the line of sight for comparison with data. 
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Table 6. Basic parameters of fits to 147 profiles of 131 GCs in NGC 5128 



Name 


Detector 


[mag] 


(V - f 606)0 
[mag] 




Model 


-> 

'I'min 


'bkg 
[Lopc^-] 


Wa 


c/n 


^0 (F606W) 
[mag arcsec"^ ] 


log ro 
[arc sec J 


log ro 
[pc] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


AATI11563 


WFC/F606 


0.308 


0.091 + 0.050 


42 


K66 

W 

S 


42.97 
42.87 
61.87 


2.56+0.15 
2.25 + 0.25 
2.57+0.17 


^ ^n+9.10 
6 10+"™ 


1 44+0.03 
'■^-0.06 
1 89+"'" 

1.92+"" 


ift #;/i+0.06 
16.64„„3 

■ 5 75+0.04 

-0.04 

15.27+"f 


1 ')nf;+0.025 
-^■206_o.oi4 
_i I2Q+0019 

^■'^^-0.020 
_2 140+0.131 


059+°"^5 
"■"^^-0.014 
135+0.019 
^■'^''-0.020 
-0 875+" '^' 


AAT113992 


WFC/F606 


0.308 


0.200 + 0.050 


26 


K66 
W 


10.19 
15.15 


-80.13+6.35 
-90.92 + 8.88 


6.60„2„ 

7 10+"'"' 
'■' -0.40 


1 4|+8:M 
'■*'-0.05 
2.46+"-33 


'*''-0.02 
16 97+" "3 
-Q.02 


-1 082+''-^'*^ 


,7oi8;4f5 
u.l 'U_Q o,4 

183+""'^ 


AATI 15339 


WFC/F606 


0.308 


0.148 + 0.050 


35 


S 

K66 

W 

S 


11.51 
47.87 
116.80 
38.38 


-73.95 + 3.30 
-42.66 + 041 
-45.30 + 2.73 
-4378 +0.92 


13 70+'-'" 

6 40+"''" 
''■*"-0.40 


1 54+8:S§ 

'■■^^-0.08 
T n8+"-23 

2 00+0.25 
2 15+^-^^ 


16.13S 
11.08+g;« 
15.78+0;'0 

1 T QT+0.26 
-0.26 


1 /;/;n+8-8'^^ 
-1.66U_oo72 
r, nio+0.191 

1 9CC+0.043 

12 511™ 
^■■"'-0.177 


-0 395+8:8W 
"■'"^-0.072 
1 7^9+0.191 
-0.214 
0|0+''"43 

'■^ -0.177 



A machine-readable version of the full Table[6]is available online (http://www.astro.keele.ac.uk/'-'dem/clusters.html) or upon request from the first author. Only a short extract from it is shown here, for 
guidance regarding its form and content. 



The complete structural and dynamical details of anvlSersid m odel 
cluster are then known in full, just as they are and 
IWilsonI models. 

Equations ^ and ([7} differ formally from equation the 
first two include a velocity scale parameter ctq, while the las t uses 
an ex plicit radial scale ro. In the formulation of his model, iKind 



1966 ) defined a radial scale associated with ctq: = 9crQ/(4nGpo), 
where po is the central mass density of the model and the numerical 
coefficients are chosen so as to make /(ro) « /o/2 in most models. 
We adopt the same definition for our single-mass, isotropic I WilsonI 
models. For both of these, we therefore have 



4nGpori 



(King and Wilson models) . 



(9) 



For the lSersij model, however, as we mentioned above it happens 
that Po ^ oo when n > I. Thus, for any n in these models we 
define a velocity scale in terms of the central surface mass density, 
Eo = T/o, and the core radius fq: 



(Sersic models only) . 



(10) 



The factor of two in the numerator here is chosen for maximum 
compatibility with equation ^ for the other models, which tend to 
have Eo ~ 2poro for the H^o or c val ues of most r eal c lusters. 

Although the model scale ro in lKind ( 1 19661) and lWilsorj mod- 
els is generally close to an observational core (projected half- 
power) radius, it is important to recognize that th e two a re not iden- 
tical in principle; an equivalence holds only for lSersij models, by 
virtue of our definition of it in equation ([8}. Similarly, the model ve- 
locity scale o"o is never equal to the velocity dispersion (projected 
or unprojected) at the centres of clusters. The connections between 
the model parameters and these observable quantities are straight- 
forward to derive for any member of our model families, but this 
can only be done numerically. 

We additionally fit all of our cluster profile s with the analyt- 
ical parametrization of I{R) introduced bv iKind ( Il962l) . and with 
"power-law" models in which I(R) oc ij'"'' at large projected radii 
and I(R) — > constant at small R. Other authors have also fitted these 
models to many old GCs and young massive clusters, but here we 
have found that doing so yields no substantiv e ly new i nfor mation 
beyond what can be learned from lKind h966h , IWilsorl and lSersid 
fits. We therefore do not discuss these other alternatives any further, 
except t o show explicitly (in 34.4| below) that they never outperform 
Iwilsorj fits to the NGC 5128 cluster profiles in any case. 



4 FITS 

The models described in 93.2l are fitted to our data after first being 
convolved with the ACS/WFC PSF for the F606W filter. Given a 
value for the scale radius ro discussed in 93.21 and some specified 
shape parameter, we compute a dimensionless model profile I^oi = 
Anod/^o and then perform the convolution 

Cd(^l'-o) = JJ T^od(R'/ro) X hsF [(x - x'), (y - /)] dx' dy' , 

(11) 

in which R^ = x- + y-; R'- = x'^ + y'^', and /psf is the PSF profile 
normalized to unit total luminosity. We use the circularly symmet- 
ric function in equation ([3} to approximate Ipsp. We also allow for a 
non-zero background /bke, and so ultimately minimize the standard 



X 



Uota(/?,)-/oX/* (i?,|ro)-/bkg 



(12) 



for the measured intensity profile and uncertainties of any cluster 
in Table|3] 

In practice, we identify the best-fit member of a model family 
by first computing unconvolved model profiles for a large n umber 
of fix ed v alues of the appropriat e shape parameter (Wo for iKind 
Il966l and Iwilsoni models; n for ISersid models). Given any one 
model in such a pre-set sequence, we convolve it with the PSF and 
then vary ro (with a new convolution required for every change), Ig, 
and /bkg until is minimized for that particular value of Wo or n. 
We do this in turn for every model in the grid to identify the single 
best fit, with the smallest from the chosen family. This proce- 
dure further allows us to estimate uncertainties for all fitted and 
derived model parameters, from the range of their values over all 
models for which is within a specified distance of the absolute 
minimum (e.g., 6x^ < 1 for 68% confidence intervals). 



4.1 Main model parameters and example fits 

Table|6]summarises the basic ingredients of all model fits to our full 
ACS cluster sample, including the duplicate profiles of the obects 
listed in Table [T] Only a sample of Table |6] is shown here; it is 
available in full, in machine-readable format, online or upon re- 
quest from the first author. 

The first column of Table |6] gives the cluster name. The sec- 
ond column reports the detector/filter combination from which we 
derived our observed density profile. This is always WFC/F606 
here, and is written out only for compatibility with analogous ta- 
bles reporting our parallel work on M31 clusters in lSarmbv et al.l 



© 2007 RAS, MNRAS 000.[Tlt26l 



Globular clusters in NGC 5128 



11 



( l2007h . The third column of the table lists the F606W-band extinc- 
tion [a constant A/rjoe = 2.8£(S - (/) = 0.308 mag for all clusters 
in NGC 5 1 28] ; the fourth column is the colour term ( V - F606)o to 
transform photometry from the native bandpass of the data to the 
standard V scale (from Table |4]l; and the fifth column records the 
number of points in the intensity profile that are flagged as OK in 
Table|3]above, and thus were used to constrain our model fits. 

The subsequent columns in Table |6] cover three lines for each 
cluster, one line for each type of model fit. Each line records: 

Column (6): identification of the model being fitted. 

Column (7): the minimum {not divided by the number of 
degrees of freedom) obtained for the best fit in that class of model. 

Column (8): the best-fit background intensity in the F606W 
bandpass. 

Column (9): the dimension less central potential Wo of the 
best-fitting model (for lKing||l966l and lWilsonl models only). 

Column (10): the concentration c = log( r,/ro) fo r lKin j ( ll966h 
an d lWilsonI models, or the index n of the best lSersi3 fit. 

Column (11): the extinction-corrected central surface bright- 
ness in the F606H' bandpass; the intrinsic K-band central surface 
brightness follows from adding the colour in column (4). 

Column (12): the logarithm of the best-fit scale radius tq in 
arcsec (see 33. 2t . 

Column (13): the logarithm of yq in units of pc (obtained from 
the angular scale assuming D = 3.8 Mpc for NGC 5128). 

To estimate the errorbars on these parameters (and those on all 
the derived quantities discussed in 34.2| below), we first rescaled the 
for all fitted models in any one family, by a common factor cho- 
sen to make the global minimum = {N^^^ - 4), where A'p,s is the 
number of points used in the model fitting. Under this re-scaling, 
the global minimum per degree of freedom is exactly one. We 
then found the minimum and maximum values of any chosen pa- 
rameter in all models with a re-scaled x^ ^ ;tmin ^- Normally, if 
were not re-scaled, this would give a 95% confidence interval on 
the parameter. Here it is just a well-defined, quantitative way to as- 
sign reasonable uncertainties to any model qu antity for any cluster. 

Figure |4] shows the best-fit iKind ( Il966h models for the GCs 
C147 and C029. The first of these is of average brightness and size 
(My = -7.6 and effective radius = 3.8 pc from the model fit). 
The second is much brighter and obviously more extended (fitted 
My = -10.2 and Ri, = 5.1 pc ). C029 was previously observed with 
STIS as part of the sample of lPaper j. where it was identified as one 
of several large clusters with extended haloes showing significant 
excess power beyond the nominal [Kinel ( 1 19661) tidal radius. This is 
also obvious in our new analysis. 

The top panels of Figure |4] show the intensity versus radius 
profiles output from ELLIPSE for each cluster, after our conver- 
sion from F6061V to extinction-corrected V, described in ^ but 
before subtracting a constant background. The lower panels of 
the figure show the intensity profiles after subtracting the fitted 
backgrounds and converting to surface brightness, jiy = IbAQl - 
2.5 log(/v/Lo pc"'). In all panels, the lower x-axis marks the pro- 
jected radius R in arcsec and the upper j-axis gives the equiva- 
lent scale in pc. The dashed curves falling steeply towards large R 
show the (arbitrarily normalized) shape of the PSF in equation l[3]l. 
The bold, d ot-dash curves are the intrinsic (unconvolved) best-fit 
lKing| ( ll966h models, with an added background in the upper pan- 
els but not in the lower panels. The solid (red) lines are the PSF- 
convolved models. Open squares are data points that were included 
in the least-squares model fitting, i.e., those flagged as OK in Ta- 
ble [3] Asterisks in the upper panels are points that were not used to 



constrain the fits (flagged as BAD, SAT, or DEP in Table |3j; they 
have been omitted altogether from the lower plots. The dotted ver- 
tical lines in all four panels are at the radius in each cluster where 
the intrinsic model intensity is equal to the fitted background level. 
For radii larger than this, the observed intensities are at least 50% 
background according to the modeling. Points with intensities be- 
low the subtracted backgrounds are represented in the lower panels 
by solid points placed on the lower x-axes, with errorbars extending 
upwards. 

The plots of the average C147 in Figure |4] are typical of the 
results for most of the clusters in our sample, in the sense that our 
derived intensity profiles extend to large enough radius that a con- 
stant background level has been reached. In this majority of cases, 
our procedure of fitting for the background at the same time as 
the intrinsic model parameters is very robust, with the same /bke 
estimated assuming any one of the three intrinsic cluster models 
discussed in 33.21 C147 is further typical of many moderate- and 
low-brightness objects with L < 10^ Lq, in that the local "sky" level 
dominates any cluster signal outside of just a few intrinsic effective 
radii. Deviations of the data from the best-fit models in any of our 
five families occur mostly in this region of noise, and thus there is 
little difference in the overall quality of fit (and no large changes in 
the derived cluster parameters) from one type of model to another. 

By contrast, the plots of C029 in Figure|4]show an interesting 
phenomenon that is seen in most of the very bright clusters i n our 
sample. As we also mentioned above, it is clear here that a lKind 
(1966,) model is simply not a good description of the data. One way 
of expressing this is to note the systematic excess of measured in- 
tensity over the best-fit model at radii > 1" 18 pc, or about 3.5 
fitted half-light radii. In order to minimize in this case, the fitted 
background intensity is spuriously high, coming in at roughly the 
average level of the outermost datapoints but failing utterly to re- 
flect the clear, continual decreasing density profile of the real clus- 
ter. A stronger way of stating this "problem" is that anv lKind i ll 9661) 
model is too concave to match the data for this cluster: the model 
does not have the shape of the observed surface-brightness profile 
even at intermediate radii, outside of the PSF-blurred core but in- 
side of where the fitted background level is reached. 

Either wa y, on the basis of this failure of the standard iKind 

( Il966l) models. iPaper 11 identified C029 and 5 other clusters similar 
to it as showing possible evidence for "extratidal light." We have 
now found several other examples like this, but we have also fit 
them with alternate structural models that have intrinsical ly more 
extended haloes by construction. We generally find that the lWilsonl 
models for such objects return lower fitted background levels, dif- 
ferent global cluster parameters such as and Ltoi, and often dras- 
tically smaller x^ value s. It would appear that the "excess" light 
implied by King j 19661) model fits to some large clusters is likely 
a symptom of generic shortcomings in the model itself — the theo- 
retical basis for it is weak in the low-density and unrelaxed farthest 
reaches of cluster envelopes — rather than the sig nature of genuine 
tidal debris. [McLaughl in & van der Marel' /2005) reached a similar 
conclusion after comparing King (1966) and Wilson fits to more 
than 150 old and young massive clusters in the Milky Way and 
some of its satellites. 

Figure|5]compares the fits of lKind jl966l) , lwilsonl and lSersi3 
models to the background- subtracted, V-band surface-brightness 
profiles of eight more globulars in NGC 5128, displayed in order 
of increasing total cluster brightness; C029 appears in the second- 
last row of this figure. The curves and points in every plot have the 
same meaning as in the lower panels of Figure |4] 

The first cluster shown in Figure[5] CI 18, is one of the faintest 
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Figure 4. Examples of I Kind il966h model fits to two clusters in our sa mple: th e fairly average C147 (left panels), and the very bright C029 (right panels), 
which has a large halo not well described by this standard model (cf. IPaper ih . Dashed curves trace the PSF intensity profile; bold dash-dot curves, the 
unconvolved best-fit model (with background added in the upper panels but not in the lower ones); and solid (red) curves, the PSF-convolved best fits. In each 
panel the lower jc-axis is marked in arcsec and the upper x-axis in parsec. Vertical dotted lines mark the radius where the intrinsic cluster intensity is equal to 
the best-fit background. See text for more details. 



in our sample {My = -6.0, or Ly ^ 2.1 x lO"* Lq according to 
the iKin^l 19661 model fit shown) and has an effective (or half-light) 
radius of =^ 3^ pc — typical of most known globular clusters in 
any galaxy — depending slightly on the model fit. The next cluster, 
C166, is similarly faint {My = -6.5, or Ly ^ 3.4 x 10* Lq) but 
significantly more extended: Rj, ^ 6.5 pc for any of the models fit. 
This conclusion is clearly not influenced by the PSF; in fact, these 
plots demonstrate that almost all of our cluster candidates are very 
well resolved indeed. 

The next two clusters, C171 and C167, both have the same 
brightness (fitted My ^ -7.4 or Ly 7.8 x lO'* Lq, essentiaOy 
at the expected peak of the GC luminosity function), but the first 
is significantly more diffuse than the second: Rj, ^ 8 pc versus 



Ri, - 3 pc. Evidently, there is substantial scatter in any size-mass 
relation that one might try to explore for average and faint globulars 
in NGC 5 128. This is reminiscent of the well-known situation in the 
Milky Way GC system, an d we retur n to the point in our analysis 
of structural correlations in lPaper ivi 

As we discussed in connection with C147 above, the three 
model fits to the globulars in this first half of Figure [5] are all very 
similar: the estimates of the sky level 1^,^^^ are all about the same, 
as are the total fitted luminosities, effective radii, and other clus- 
ter parameters. The minimum for the different models are also 
very si milar fo r any one cluster, and there is no sign that the ex- 
tended [^Iso^ or lSersid haloes are sy s temat ically either preferred 
over or bettered by the standard [Kind i ll 9661) structure. It is worth 
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Figure 5. Examples of lKind jl966l) . lwilsoi]l and lSersi3 model fits to 8 globular clusters in NGC 5128. Results are shown in the same background-subtracted 
format as the lower panels of Figure |4] with all points and curves having the same meaning here as there. Clustercentric distance is marked in arcsec along 
the lower .v-axis in each panel, and in parsec along the upper jc-axis. Clusters are presented in order of increasing brightness from top to bottom: My = -6.0 
for CI 18, Mv = -6.5 for C166, My = -7 .4for C171 Mv = -lA for C167, My = -8.7 for CI 13, My = -8.8 for F1GC15, My = -10.2 for C029, and 
My = -11.2 for C007 (all magnitudes from lKin J 19661 model fits and assuming D = 3.8 Mpc). 



© 2007 RAS, MNRAS 0Q0.[Tll26l 



14 D. E. McLaughlin et al. 



1 10 100 



R/pc 

1 10 100 



1 10 100 




0.1 1 10 0.1 1 10 0.1 1 10 

R/ arcsec 



Figure 5 - continued Examples o flKingl h966l). IWilsorl andlSersid model fits to 8 globular clusters in NGC 5128. Results are shown in the same background- 
subtracted format as the lower panels of Figure |4] with all points and curves having the same meaning here as there. Clustercentric distance is marked in 
arcsec along the lower A:-axis in each panel, and in parsec along the upper .v-axis. Clusters are presented in order of increasing brightness from top to bottom: 
My = -6.0 for C118, My = -6.5 for C166, My = - 7.4forC17 1, My = -7.4 for C167, My = -8.7 for C113, My = -8.8 for F1GC15, My = -10.2 for 
C029, and My = -11.2 for C007 (all magnitudes from lKind[T963 model fits and assuming D = 3.8 Mpc). 
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Table 7. Derived structural and photometric parameters from 147 profiles of 131 GCs in NGC 5128 



(1) 



Detector Model log log Rc log Rji 

[pc] [pc] [pc] 

(2) (3) (4) (5) (6) 



log(R,,/R,.) 



log /o 
(8) 



log ia 
(9) 



log Ly 

[io.v] 
(10) 



[mag] 
(11) 



log '* 

(12) 



AATl 11563 


WFC/F606 


K66 


1 SO+O'" 


n35+""22 
"■"^^-0.012 


n dfiS+0.005 
0.468„j„„ 


Q 433+0.017 
"■*^^-0.032 


r, 07+0.02 
^■°'-0.03 


3 53+0.03 

3 42+8:83 
^■^^-0.04 


5 05+" "2 
^■"^-0.02 


20 10+""' 
-Q.Q5 


3 32+0.02 
^■^^-0.02 






W 


2 03+""* 
-0.08 


QOgn+0.015 
"■""'-0.016 


486+""" 
"■^°"-0.014 


3q6+" "33 
"■^'"-0.029 


3 07+0.03 
■ -0.03 


5 07+" "2 
^■"'-0.02 


20 05+""6 
■^-0.06 


3 30+0.03 
^■^"-0.03 






S 




-0 875+" '^' 


Q 478+0.010 
"■^"^-0.009 


1 353+0.155 
'■^^^-0.140 


4 42+" "9 
^■^^-0.08 


4 23+"^^' 
^■^^-0.19 


5 05+" "2 
^■"^-0.02 


20 10+""' 
-Q.Q5 
20 02+""' 
^"■"''-0.09 


3 30+" "2 
^■^"-0.02 


AATl 13992 


WFC/F606 


K66 


1 ^8+0.07 


n |44+8:4(ti 

-0.011 


560+""'" 
"■■^ ""-0.023 


Q 416+0.051 
"■^'"-0.030 


3 70+" "2 
-0.02 


3 75+0.03 
■ -^-0.02 


5 ns+0.04 
5.08„j,3 


3 l6+" "3 
^■'"-0.05 






W 


2 64+0.31 


151+""12 
"■"'-0.011 


749+"''" 
"■ '^'-0.097 


Q 559+0.155 
"■^''-0.109 


3 69+" "2 
^■"'-0.02 


3 24+" "3 
■ -0.03 


5 21+0.08 


19 69+""' 
"■"'-0.21 


2 92+"13 
^■'^-0.21 






S 




-0.395+""'' 


n cTQ+0.016 
0.528„j„5 


923+0.088 
"■'^^-0.086 


4 03+""" 
^■"^-0.04 


3 45+0.10 
■ -^-0.09 


5 05^^ 
■^■"■^-0.02 


20 11+""'' 
^"■' -0.06 


3 19+0.03 
■ -0.03 


AATl 15339 


WFC/F606 


K66 


1 ^0+0.02 
'■^^-0.01 


1 7r',+8:?^o 

'■"■'-0.214 


Q 376+0.003 
"■^'"-0.000 


T 127+0.216 
''■''' -0.190 


6 07+"^™ 
"■"'-0.17 


7 52+"^"^ 


5 1 9+0.02 
■^■"-0.02 


19 76+0.05 
-Q.Q5 
19 67+" "7 
'■ -0.08 


3 64+"^"^ 
^■°*-0.02 






W 


2 01+0.20 


-0 033+""35 
"■"^^-0.040 


397+""44 


430+""*" 
"■'^"-0.062 


4 1 9+0.05 
'■' -0.04 


3 97+8:55 

-0.08 


5 27+0.03 
■ -0.03 


3 63+"-"" 
^■"^-0.06 






S 




-1 246+" 


384+'i-''^5 
"■^°^-0.009 


1 630+"- '86 
'■"^"-0.182 


4 93+0.11 
^■'^-0.11 


5.08!B;| 


5 20+" "2 
-■■^"-0.02 


19 72+0.05 
"■' -0.05 


3 64+" "2 
^■" -0.02 



A machine-readable version of the full Table[7]i.s available online (http://www.astro.keele.ac.uk/--'deni/clusters.htnil) or upon request from the first author. Only a short extract from it is shown 
here, for guidance regarding its form and content. 



Table 8. Derived dynamical parameters from 147 profiles of 131 GCs in NGC 5128 



Name 


Detector 


^pop 


Model 


log Mioi 


log Ej, 


log So 


logpo 


log £(, 


log CTp 


log I'csc.O 


log f,t 


log h 










[Mo] 


[erg] 


[Mopc-2] 


[Mg pC-3] 


[Me pc-2] 


[kms-'] 


[kms-'] 


[yr] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


AATl 11563 


WFC/F606 


1 939+0.239 
'■"'-O240 


K66 

W 

S 


5 34+0.05 
^■■'*-0.06 
5 36+" "' 
■^■^"-0.06 
5 34+0.05 


=.f\ /:7+0.08 

50.62„„, 

50 6''+""* 
■^"■"--0.09 

50.62+""* 


4 16+""' 
^■'"-0.07 
4 1 1 +0.06 
'■"-0.06 
470+0.10 


1 07+0.06 

■^■°''-0.08 
3 7 1 +0.06 
^■"-0.07 
452+0.21 

3 71+8:4^ 
3 70+8:8? 

^■'"-0.07 
391+0.10 


1 fin+0.06 
^■S''-0.06 
3 59+0.06 
■ -0.07 
3 58+0.06 


n 9^6+0.027 
0.836„„3|, 

839+" "" 

■ -0.030 
Q 797+0.028 

0893+8:853 

895+'^-^^^ 
"■^^■^-0.036 
Q 877+0.032 

1 047+''-^^^ 
'■""-0.029 
0972+""26 

913J8:83? 

■ -0.031 


1 421+0.027 

■ -0.031 
1 427+0.027 

■ -0.031 
1 470+0.028 


8 92+0.03 
''■^^-0.03 
8 96+0.04 
"■^"-0.04 
8 94+0.03 


079+""'* 
"■"'^-0.056 

-0 049+"""^ 
"■"'^-0.045 

1 134+8.287 

_0 i97+8:S35 

_0 215+8:819 
"■'"■'-0.046 
158+""' 

3 482+^-'^^^ 
^■'"^-0.391 

067+""** 


AATl 13992 


WFC/F606 


2 914+0.431 
^■"+-0.423 


K66 

W 

S 


555;8:85 

5 68+8:?^ 

5 51+0.06 


5Q 94+8:8? 

51 01+'^-"^ 
5Q 92+0.09 


4 16+8:48 

4 16+'^-'^^ 
'■'"-0.07 
449+0.07 


363^8:8^ 

3-38+Q^27 
3 66+""' 
3 94+0.05 
^■'^-0.06 
3 94+0.06 


1 476+8:85i 

1 499+8:81 
'■'"-0.036 
1 507+0.032 

1 767+8:85§ 
'■'"'-0.030 
1 563+""^* 


9 16+8:8? 

9 50+8--^ 
'■■^"-0.17 

9 10+""" 
■ -Q.Q5 


AATl 15339 


WFC/F606 


2 015+O238 
■ -0.235 


K66 
W 


5 49+8:8? 
■^■^'-0.06 
5 53+0.06 


51 ns+8:i^ 
51.08„„g 

51 05+""' 
■ -8BS 


6 37+8:§5 
"■^'-0.18 
4 49+0.07 


7 07+8:15 

4.22+"'" 


8 oc+0.03 
8 89+0.08 








s 


5 5i;8:8§ 

■^■■^'-0.06 


51 05+""° 
■"■"■^-0.08 


2:8:!? 
■ -0.12 


c 78+8:5^ 

^■^'*-0.25 


3 94+8:8? 
^■'^-0.06 


■ ;8:859 

'■"' "-0.030 


8 07+8:8^ 
'^■*"-0.03 


°- I:8:S§§ 

■ -0.343 



A machine-readable version of the full Tablefslis available online (htlp://www.astro.keele. ac.uk/-dem/clusters. html) or upon request from the first author. Only a short extract from it is shown here, for 
guidance regarding its form and content. 



noting that the ISersid model fits can return significantly smaller 
core radii and brighter intrinsic central surface brightnesses than 
the other models; but at some level this simply points to a limita- 
tion of the data, since the unavoidable PSF convolutio n eflrec tively 
erases the intrinsically cuspy inner structures of n > l lSersid mod- 
els, which might be disallowed by much higher-resolution data. 

In the second half of Figure[5] the clusters C113 and F1GC15 
are both taken from the bright side of the GC luminosity function 
in NGC 5128 (My = -8.7 and My = -8.8), and once again one 
is much larger than the other {R/, = 2.5 pc against Ri, = 16.5 pc in 
the most conservative, lKing|[l96^ model fits). Now, however, the 
quality of fit difi'ers significantly between the different mo dels. In 
CI 13, the m inimum i s nearly 7 times smaller for the IWilsonI 
model versus Kine ('1966*), and just over 5 times smaller for lSersij 
versuslging ( 1966). F orFlG C15, the Wilson fit has a x' about 20% 
larger than the lKingl jl966h fit (so the latter is formally somewhat 
better, but only marginally so), while the iSersic model is effec tively 
ruled o ut with a minimum more than 6 times larger than the lKind 
( Il966l) model. 

The last two clusters illustrated here are C029 — confi rming 
that th is object is much better fit by a lWilsonI model than a lKind 
(119661) model, and also showing that the former does better than a 
ISersid model — and C007, which is the brightest GC in our sample 
(My = - 1 1.2, or L y = 2.6 x 10* Lq from a lKing 1966 model fit). It 
too is better fit by a IWilsonI model than by either of the other two. 

We note again that, because the extended-halo structure of all 
the bright cluste rs in Figure |5] is poorly fit by the relatively com- 
pact lKinj d 19661) models, the fitted background intensity is notice- 
ably, but spuriously, higher in the leftmost panels than in the other 
columns. 



We will return to the question of which model best fits the 
majority of GCs in NGC 5128, in 94.41 below. First, we present a 
number of physical cluster properties derived within the framework 
of each of the model fits summarized in Table[6] 



4.2 Derived quantities 

Table [7] contains a number of other cluster properties derived from 
the basic fit parameters given in Table [6] Columns following the 
GC name, detector/filter combination, and fitted model are: 

Column (4): log r, = c + log rg, the model tidal radius in pc. 



which is always infinite for lSersiq models. 

Column (5): log R^, the projected core radius of the model 
fitting a cluster, in units of pc. This is defined by I(Rc) = /o/2, or 
li(Rc) - fiQ + 0.753 and is n ot nece ssarily the same as the radial 
scale ro in Table[6l except for lSersij models. 

Column (6): logR/,, the projected half-light, or effective, ra- 
dius of a model. Half the total cluster luminosity is projected within 
Ri,. It is related to rg by one-to-one functions of Wq (or c) or n and 
is reported here in units of pc. 

Column (7): log(Ri,/Rc), a measure of cluster concentration 
that is somewhat more model-independent than c or n, in that it 
is generically well defined in observational terms and its physical 
meaning is always the same (cf. our earlier discussion in 93.2b . We 
consider it a more suitable quantity to use when intercomparing the 
overall properties of clusters that may not all be fit by the same kind 
of model. 

Column (8): log Ig, the best-fit central (R = 0) luminosity sur- 
face density in the V band, in units of Lo.vpc"-. This is obtained 
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Table 9. Galactocentric radii and ^■-space parameters from 147 profiles of 131 
GCsinNGC5128 



Name Detector 
(1) (2) 



[kpc] 
(3) 



Model 
(4) 



(S) 



(6) 



(7) 



AATl 11563 WFC/F606 11.93 K66 
W 
S 

AATl 13992 WFC/F606 4.26 K66 
W 
S 

AATl 15339 WFC/F606 3.82 K66 
W 

S 



-0.608_„„^ 
_n TO I +0.040 
n ft^-7+0.039 

-0 7,26*^'^^^ 
"■^^"-0.084 
_05Q7+0.045 

-0 375+8:853 
"■^"-0.041 
-0.467+°"™ 

""■^=*'-0.042 



4.660;"»* 
4 642+8:8?5 
+■"^^-0.079 
A /;m+0.070 
4.607_„j„g 

4559+0.084 

44in+8:l^i 
'"-0.246 

4 714+0.081 

= 147+8:8?^ 
-0.071 
5.070;"'"5 

5 (139+8:8^2 
■^■"^^-0.074 



n 34Q+0.004 

n ^n(;+0.008 
0.306 „ ™, 

O345+O.0B5 
n 380+8:8^5 
0.328+9™ 



'-0.012 



A machine-readable version of the full Table [9] is available online 
(http://www.astro.keele.ac.uk/-denVclusters.html) or upon request from the first author. Only a 
short extract from it is shown here, for guidance regarding its form and content. 



from the fitted central surface brightness in Column (11) of Ta- 
ble|6] by first applying the "V-colour" correction in Column (4) of 
that table to obtain the central jiyQ, and then using the definition 
log ly = 0.4(26.402 - yUv), where the zeropoint corresponds to a 
solar absolute magnitude of Mq.v = +4.83. 

Column (9): log /'p, the V-ba nd l uminosity volume density 
(Lo.v pc"^) at r = for lKind ( 1 19661) andlWilson ' models but at the 
three-dimensional radius r - ro for lSersij models. In the first two 
cases, id = Srio/ro where J7" is a smooth, model-dependent func- 
tion of Wq or c, which we have calculated numerically. For lSersi3 
models, as we discussed in 33.21 the unprojected density is infinite 
as r ^ when n > I, and thus we only quote jo = /Co) for any 
fitted n. This finite quantity is related to the quotient lo/rg by a well- 
defined function of n, which we have again computed numerically. 

Column (10): log Ly, the total integrated model luminosity in 
the V band. It is related to the product lor^ by model-dependent 
functions of Wo or c, or n. 

Column (11): V,„t = 4.83 - 2.5 log (Lv/Lo) + 5 log(O/10pc), 
the total, extinction-corrected apparent V-band magnitude of a 
model cluster, assuming D = 3.8 Mpc. 

Column (12): log/;, = log(Lv/2nRl), the V-band luminosity 
surface density averaged over the half-light or effective radius, in 
units of Lo.vpc"". The V surface brightness averaged over R/, is 
(^iv)i, = 26.402- 2.5 log (/a/Lo.v pc-^). 

The uncertainties on all of these derived parameters have been 
estimated from the re-scaling procedure described after Table 
|6] If the distance D to NGC 5128 is difi'erent from our adopted 3.8 
Mpc, the quantities in Table|7]will change according to: rm, R^, and 
Ri, all oc D; Iq and /;, independent of D; jo f ' ; and Ly oc D-. 

Table [8] next lists a number of cluster properties derived from 
the structural parameters already given plus a mass-to-light ratio. 
The first two columns of this table contain the name of each cluster 
and the combination of detector/filter for our observations of it, as 
usual. Column (3) lists the K-band mass-to-light ratio that we have 
adopted for each object from the analysis in 33.11 assuming for 
definiteness an age of 13 Gyr for all clusters. The errorbars on T^"'' 
in Table [8] are larger than those in the 13-Gyr column of Table [5] 
as we allow now for a ±2-Gyr uncertainty in age on top of the 
previously tabulated uncertainties in [Fe/H]. The remaining entries 
in Table[8]are, for the best fit of each model to every cluster: 

Column (5): log Mto, = log T^"'' + log Ly, the integrated 
model mass in solar units, with log Ly taken from Column (10) 
of Table|7] 

Column (6): log Et, the integrated binding energy in ergs. 



defined through Ei, = -(1/2) J^' Anr^p(fidr. Here the minus sign 
makes Ei, positive for gravitationally bound objects, and 0(r) is the 
potential generated (through Poisson's equation) by the model mass 
distribution p(r). E), can be written in terms of the fitted central lu- 
minosity density jo; scale radius ro; a model-dependent function 
of Wo o r c, or n; an d Y^"''. A more detailed outline of this proce - 
dure for lKind (Il966i) models may be found in iMcLaughlinI ( |2000|) . 
which we have followed closely to evaluate Ei, for our other model 
fits as well. 

Column (7): log Eq = log T^"'' -I- log /q, the central surface 
mass density in the model, in units of Mq pc"^. 

Column (8): log po = log T^"'' + log /o , the central volume 
density in units of MqPC"^ (except for lSersig models, where po is 
the density at the three-dimensional radius ro). 

Column (9): log E;, = log T'^"'' + log /;,, the surface mass den- 
sity averaged over the inner effective radius Rj,, in units of Mq pc"-. 

Column (10): log iXpo, the predicted line-of-sight velocity dis- 
persion at the clu ster cen tre for .Kina t l966) and Wilson models, 
but at ffl = Rc for lSersij models, in km s The solution of Pois- 
son's and Jeans' equations for any model yields a dimensionless 
(Tp o/cTo, and with (To given by the fitted ro and po = T^°''7o through 
equation ([9]l or l IlOK the predicted observable dispersion follows 
immediately. 

Column (11): log Ves^.o, the predicted central "escape" velocity 
in km s"' . A star moving out from the centre of a cluster with speed 
VescO will just come to rest at infinity. In general, then, 

Vesc.o/c^o = 2[Wo + GM,„,/r,tT2] . 

A (finite) dimensionless Wo is associated with any n > for lSersid 
models by solving Poisson's equation with — > in the limit r — > 
oo. The seco nd term on the right-hand side of the definition of Vj-s^.o 
vanishes for lSersi j models, in which r, — > oo. 

Column (12): log ?,h, the two-body relaxation time at the 
model projected half-mass radius, in years. This is estimated as 

_ 2.06 X 10^ yr M,'„f /;,f 
* ln(0.4Mto,/)!u) m^, 

teinnev & Tremain3ll987l . eq. 8-72), if m-i, (the average stellar 
mass in a cluster) and Mmt are both in solar units and R/, is in pc. We 
have evaluated this timescale assuming an average m^, = 0.5 Mq in 
all clusters. 

Column (13): log/o = log Jpo/(27rt7"^)''^j, a measure 
of the model's "central" phase-space density in units of 
Mq pc"^^ (km s"')"'. In this expression, cr^ refers to the central one- 
dimensional velocity dispersion not projected along the line of 



© 2007 RAS, MNRAS 000.[TH26l 



Globular clusters in NGC 5128 



17 



sight. The ratio (tJctq is obtained for iKind ( Il966h and Iwilsoni 
models from the solution of the Poisson and Jeans equations for 
given W(i or c, and ct q is kii own from equation ([9}. This proce- 
dure breaks down for ISersid models in general, since the unpro- 
jected p — » oo and cr —> as r —> for n > 1. As before, then, 
for these models we instead calculate /o at the three-dimensional 
radius r - rQ. In any case, with the centra l relaxation time t,c 
of a c luster defined as in equation (8-71) of Isinnev & Tremaind 
( Il987h . taking an average stellar mass of m^, = 0.5 Mq and a typi- 
cal Coulomb logarithm In A a; 12 leads to the approximate relation 
log fe/yr)- 8.28 -log /o. 

The uncertainties in these quantities are estimated from their 
variations around the minimum of on the model grids we fit, as 
above, but now combined in quadrature with the uncertainties in 
the population-synthesis model Y^"''. Note that if T^°'' is changed 
to any other value for any cluster, or if any distance other than D = 
3.8 Mpc is adopted for NGC 5128, the properties in Table[8] scale 
as M,ot oc TZ)2; Eb oc Y- D^; So TZ)" and E,, oc T/)"; po oc 
Y£»-'; o-po oc T'/2z)i/2 and Ve,co oc T"- D"-; t,t, oc T'/2£)5/2. an^j 
/oocT-i/2z)-5/2. 

Table|9]provides the last few parameters required to construct 
the fundamental plane of globular clusters in NGC 5 128, under any 
of the equivalent formulations of it in the literature. The first of 
these remaining parameters is the projected galactocentric distance 
ifgc, in kpc, of each cluster. This is listed in Column (3) of Table 
[21 following the cluster name and detector/filter combina tion. It is 
obtained from the angular R,^ given in Tables 1 and 2 o f IPaper ij. 
assuming D = 3.8 Mpc as always. 

The last three columns of Table [9] give analogues of the 
so-ca lled k parameters introduced by Bende r. Burstein. & Fabed 
( Il992 ). who define an orthonormal "coordinate" system based on 
the three main luminosity observables of galaxies. Similarly to 
them, but in terms of mass-based quantities instead, we define 

KmA = (log Q + log Rl,)/ V2 

K„,2 = (log 0-2 ^ + 2 log E,, - log Rl,)/ V6 (13) 
a:,„.3 = (log (T^ - log E/, - log R,,)/ V3 

which differ from the k's in ISender. Burstein. & Faberl (and 
iBurstein et al.lll997h only in our use of the average mass surface 
density X;, rather than the luminosity surface density /;, = E;,/Tv. 
The main reason for doing this is to facilitate the comparison of 
fundamental planes for young and old star clusters, by removing 
the influence of age on /;, for two objects of similar mass and size. 
It similarly avoids the influence that any metallicity dependence in 
Yv has on a luminosity-based fundamental plane. 

In equation l |13l l, k„,j «-> log (cr^gRi,) is related to the total 
mass of a system, and k„,,3 <-> log((r^gRi,/Mf„t) contains the ex- 
act details of this relationship — that is, information on the inter- 
nal density profile p( r). In f act, the mas s-based k„^3 can be viewed 
as a re placement for iKiiig!- or IWilsoiil -model concentrations c or 
ISersid -model indices n, or any other model-specific shape param- 
eter. As such, any trends involving k^^t, are directly of relevance 
to questions concerning cluster (non)homology. The definition of 
Km.2 <-» log (LI) is chosen to make the three Km axes space mutually 
orthogonal in the parameter space of ((t^q,Ri,,'I,i,). 

In calculating k„, i, k„,2, and k„,j for Table|9] we have used the 
(Tp o predicted in Table[8](Column 7) by our adoption of population- 
synthesis mass-to-light ratios. Similarly, E;, is taken from Column 
(9) of Table [8] As a result, our values for /c,,, 3 are independent of 
T^"''. Rl: is taken from Table [7] but put in units of kpc rather than 




10= 10' 20 18 16 14 13 10 

Lj/Lg [central] /i, „ [central] 

Figure 6. Comparison of concentration parameters, intrinsic radii, total 
luminosities, and intrinsic central surface brightnesses for fits of isotropic 
Iwilsoni Il975l) models to each of the two independent intensity profiles mea- 
sured for 15 of the GC candidates listed in Table[T](not including the pos- 
sible star, C156). In the upper-right panel, filled symbols refer to estimates 
of the model scale radii ro for the clusters, while open circles refer to the 
eft'ective (projected half-light) radii Ri,. 



pc, for be tter compatibility with the origina l, galaxy-oriented defi- 
nitions of lBender. Burstein. & Fabeil l ll992l) . 



4.3 Consistency checks 

There are two ways in which we can assess the internal consistency 
of the various model fits we have performed. 

First, Figur e [6] compares some of the main parameters ob- 
tained from Wilson fits to the 15 clusters in Table [T] for which we 
have two independent surface-brightness profiles measured on two 
different main ACS fields. (We have excluded the object CI 56 from 
this comparison; cf. Table|2]) 

The X-axis of every panel in Figure |6] marks the parameter 
value measured for the clusters in the survey target fields where 
they lie closest to the chip centre (the fields highlighted in bold in 
Table [T). The y-axes then measure the parameter values found in 
the ACS field where the clusters lie further from the centre. There 
is evidently no substantial or systemat ic scatte r around the line of 
equality drawn in each case. For these Iwilsoni fits we find that the 
average and rms scatter in the difference of fitted concentrations is 
(in the sense [off-centre] minus [central]) (Sc) = -0.02 ± 0.08, as 
compared to an rms errorbar of ± 0.13 on the individual c values. 
For the scale and half-light radii, (6(\og ro)) = {6{\og Rj,)) = 0.01 ± 
0.03 dex, versus an rms fit uncertainty of about ± 0.04 dex in both 
cases. Finally, we have (i5(log Ly)) = ± 0.01 dex and (Siiv.o) = 
0.03 ±0.11 mag arcsec"-, against rms errorbars of ±0.03 dex and 
± 0.09 mag arcsec"^. Nearly identical results hold for comparisons 
of fitted King ( 1966) and Sersic model parameters. 

We conclude that any PSF variations over the ACS/WFC chip 
have not introduced any systematic errors into our model fitting, 
which assumed a single PSF shape in all cases. For all subsequent 
analyses, we include the clusters with double measurements by 
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Figure 7. Comparison of total lKin3 j 19661) and lwilsot] model V magni- 
tud es versus aperture or PSF estim ates (extinction-corrected) from lPaper ij 
(cf. lHarris. Harris. & Geislej2004 . for 124 reliable GC candidates. 



adopting the fit parameters obtained from the image on which the 
cluster is closest to the chip centre. 

Second, Figure|7]shows the total apparent V-band magnitudes 
inferred from our King ( 1966) and Wilson model fits (from Col- 
umn [11] of Table |4j to the V-band equivalents of Washington 
Ti magnitudes estimated from gro und-based aperture or PSF pho- 
tometry and tabulated in [ Paper III (numbers taken originally from 
iHarris. Harris. & Geisleill2004l) . To transform Ti to V, we first cor- 
rected for extinction acc ording to An = 2.54 x E{B - V) = 
0.279 mag for all clusters jHarris &Cantemall 19791) . and then used 
the relati onship Koi = J i.o + 0.052 -I- 0.256 x (C - TOo from lPaper j 
(s ee alsolGeisleill 19961) given the de-reddened Washington colours 
in IPaper iJ and our Table |4| above. We have excluded from this 
comparison the clusters C168, F1GC34, and F2GC14, for the rea- 
sons given in Table |2l the objects C145, C152, C156, which are 
likely foreground stars; and C177, which is probably a background 
galaxy. 

For bright clusters with Vto( < 20 (My < -7.9), our model 
luminosities compare rather well in general to the aperture or PSF 
magnitudes; with (5Vtot = [Vtot(tnodel) - Vtot (aperture o r PSF)], we 
have <(5V,„,) =0 + 0.13 mag (rms) for the 'King' h96&) model fits, 
and {SVtat) = -0.14 ±0.18 mag (rms) for Wilson fits. Both results 
are quite acceptable, considering (1) the number of steps involved 
in transf ormi ng either our models or th e Wash ington photometry of 
IPaper III and iHarris. Harris. & Geisleil ( |2004|) to standard V mag- 
nitudes; and (2) the fact that the aperture or PSF magnitudes are 
from analyses of ground-based data and limited to clustercentric 
radii of only R < 2"-3" (about 35-55 pc), which can easily miss 
up to ~ 0.2 mag of cluster ligh t (see Figure [5] above, and also the 
discussion in lHarris. Harris. & Geisler 2004). The 0.15-mag offset 
between the average 'Wilson and King ( 1966) model magnitudes is 
a direct result of the fact that the latter models are generally too 
compact to fit the brightest clusters well and thus underestimate the 
true luminosity of the cluster haloes. 

Fainter clusters with > 20 according to the aperture or 
PSF measurements have <5Vtot olfsets that are on average ^ 0. 1 mag 



brighter than for the brighter GCs, and that have a larger rms scat- 
ter about the average. This is due to our ability here (unlike in 
IHarris, Harris. & Geislerl [20041) to resolve all of these faint clus- 
ters and robustly to estimate and correct for the local background 
intensity around each of them. We thus view our present measure- 
ments of Vtot from model fits to resolved surface-brightness profiles 
as being more reliable than the previously published numbers. 

A comparison between our Sersic-model Vtoi and the ground- 
based estimates gives re sults that are nearly indistinguishable from 
those shown here for the lKind j 19661) fits. 



4.4 Quality of fit for different models 



We now return to the issue, raised in 34.11 of deciding which model 
family bes t describes the structure of GCs in NG C 5128 in general. 
FoUowing lMcLaughlin & van der Marell ( |200^ . we define a statis- 
tic that compares the x' of the best fit of an "alternate" model for 
any object to the;^'^ of the best-fit standard King (1966) model: 



' C\f alternate ^K66)/C\'alternate XvidS^ ■ 



(14) 



This statistic vanishes if the two models fit the same cluster equally 
well. It is positi ve (with a m aximum of -I- 1) if the alternate model is 
a worse fit thanl King| ( IT966h . and negative (with a minimum of -1) 
if the alternate model is a better fit. We have found that for fairly 
modest values, -0.2 < A < 0.2, even though one model formally 
fits better than the other the dilference is generally not significant. 
For A more towards the extremes of ±1, the improvement afforded 

by the model with lower is more substa ntial. 

Figure [8| shows this A statistic for Wi lsonI - and ISersi^ model 
fits (open square s and asterisks, respectively, in all panels) ver- 
sus iKingj h966l) fits to 124 GCs in our current sample. These 
numbers are plotted against the clusters' half-light radii Ri„ total 
model luminosities Ly, and the intrinsic average surface brightness 
{jiv)h = 26.402 - 2.5 \og{Lvl2nRl) as estimated fr om the Kiii^ 
( Il966l) fits. 

The left-hand panel of Figure[8|shows immediately that many 
cl usters in N GC 512 8 are better described by th e distended shapes 
of lWilsonl or lSersid models than by the concave iKind ( 1 1966 ) pro- 
files. However, the intrinsic size of a GC is clearly not a good pre- 
dictor of which model shape will fit it best. 

The middle plot, of A vs. Ly, shows that clusters brighter than 
Ly ~ (1-2) X 10' Lq — that is, from just above nominal turnover 
point of the GC luminosity function — are generally fit much better 
by models with more extended haloes than allowed in the stan- 
dard lKing| ( IT96^ descrip tion. Fainter clusters are u s ually fi t almost 
equally well by WilsonI models as bv iKind i ll966l ). [Sersid models 
tend to fare si milarly, alt hough there are a number of individual 
cases in which lWilsonI and lKind ( Il966l) models both do better. 

The graph of A versus {iiy)i„ in the right-hand panel of Fig- 
ure [8] is the most telling. This shows clearly that clusters with 
high average surface brightnesses are significantly better fit by 
models with larger haloes than King ( 1966). At some level this 
might seem almost a tautology, given that A correlates with Ly 
and that {jiy)], obviously does as well. But it is interesting at 
this point to compare with the si t uation for Local Group clus- 
ters. iMcLaughlin & van der Marell ( |2005|) show that 153 globu- 
lars and young massive clusters in the Milky Way, the LMC and 
SMC, and the Forn ax dwar f spheroidal are systematical ly better 
fit by the extended IWilsorj haloes than by iKind l ll966l) models 
whenever the clusters' surface-brightness profiles are accurately 
defined to large projected radius (more than about 5 Rh) such that 
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Figure 8. Relative quality-of-fit, A (eq. 1141 ). for isotropic lWilsotJ and lSersi j models (open squares and asterisks) versus iKindiiged) models, for 124 GCs in 
NGC 5128. (Objects C168, F1GC34, F2GC14, C1 45, C15 2, C156, and C177 are excluded from the starting sample of 131 cluster candi dates for the reasons 
given in Table|2]) The right-hand panel shows that lwilsoni models, with their relatively more extended haloes, tend to give better fits than lKindfT96^ models 
for high surface-brightness GCs, which have envelopes that are empirically better defined against the background sky out to large clustercentric radii. 




Figu re 9. Lef t: Rel ative q uality-of-fit statistic A (eq. ITU) for lSersid versus lKin3 h96(j) models fitted to 124 rehable GC candidates in NGC 5128, against 
A for I Wilson v ersus fits to the same clusters. Right: Analogous comparison for power-law versus [wilson models, against that for King ( 1963) 

versus I WilsonI models. L abels and bold lines indicate which of the three models being compared is the formal best fit for GCs falling in different regions of 
either graph. An isotropic Iwilsod j 19751) model nearly always fits at least as well as any of the others we have tried, and often does significantly better. 



the intrinsic structural differences between the model haloes be- 
co me rehably observable in the first place. T here is no corr elation 
in iMcLaughlin & van der MmbII between the IWilsonl - fKind A and 
cluster luminosity or mass. In our case, clusters with brighter {^v)h 
are higher above the background in NGC 5128 over a larger frac- 
tion of their volume, and thus their far halo regions are necessarily 
better defined in purely ob servatio nal terms. It is when this happens 
that th e large envelopes of lWilsor] models particularly (and those of 
ISersid models to a lesser extent) are preferred over the more com- 
pact iKind (Il96^ structure. Fainter and lower surface-brightness 
clusters may very well be similarly extended in general, but the 



data at large radius in these cases are not good enough to show this 
definitively. 

We therefore conclude from Figure [8] that the haloes o f GCs 
in NG C 5128 are generically more extended than the classic iKind 
( Il966l) model allows. These globulars thus appear to be very much 
like those in the Mi lky Way and some of its satellite galaxies (see 
iBarmbv et al ] |2007l for a discussion of GCs in M31, where the sit- 
uation is less clear). Any correlation between A and total cluster 
luminosity appears to derive from this more basic fact; it does not 
imply that massive clusters have extended haloes just because they 
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are massive per se, nor that bright clusters have systematically dif- 
ferent intrinsic structures than faint ones. 

The slight collective offset of the asterisks upwards from the 
squares in Figure [8] suggests that the glob ulars in NGC 5 1 28 are 
generally somewhat better fit by Iwilsorj models than by ISersid 
models. This is confirmed in the left -hand pa nel of Figure|9l which 
shows the individual A statistics for lWilsonl -model fi ts to th e same 
sample of 124 GCs as in Figure[8l against the A for lSersid -model 
fits. Points in the upper-right hand quadra nt (Aw > a nd As > 0) 
of this plot represent clusters for which a lKind jl966h model fits 
better than either of the other two. Clusters that are formally fit best 
by a Sersic model have As < and As < Aw, and thus fall within 
the horizontal half-trapezoidal region (defined by bold lines) in the 
lower part of the plot. When a .Wilsort model fits best (Aw < and 
Aw < As), the clusters populate the vertical half-trapezoidal area 
on the left-hand side. 

About 60% (75/124) of our clusters fall in this latter region 
of the l eft-hand panel in Figure |9] and are thus formally best fit by 
IWilsod models. Fully 90% (111/124) fall within the shghtly more 
gener ous limit s Aw < 0.2 and Aw < As -(- 0.2, which include cases 
where Wilson models are just not significantly worse than either 
[Sersic or King (1966) models. Presumably these results reflect that 
GCs have relaxed, n early i sothe rmal core struct ures tha t are cap- 
tured by [Wilson andKin^ ( 1 1966) models but not lSersi3 models in 
general , plus s trong halo components t hat are better matched by 
IWilsonI or lSersi c models than by those of King ("l96f)). 

Points well inside the horizontal half-trapezoid in the lower 
part of the left-hand panel of Figure [9] represent clusters that are 
noticeab ly better fit bv lSersij models than by either lKing| ( [T96^ or 
IWilsonI models. Although there are relatively few such objects — 
only 12/124 have both As < -0.1 and (As - Aw) < -0.1 — it might 
be thought that they could be of particular interest. If, for example, 
they were all extremely massive, this might then be supportive of 
suggestions that the brightest globulars were stripped galaxy nu- 
clei or possibly ultra-compact dwarfs. However, the average total 
luminosity of the 12 clusters just mentioned is a completely nor- 
mal (Lv) - 5 X 10"* Lq — nearly at the peak of the GC luminosity 
function — and all have Ly < 2 x 10' Lq. 

Finally, the right-hand panel of Figure |9] shows the relative 
quality of fit for two other models that we fit to all of our ACS 
cluster data: the analytical l(K) models o f lKingHl962l) , and power 
laws with constant-density cores (see the end of 33. 2t . Now, how- 
ever, we have formed A statis tics by comparing the best-fit o f 
these models to those of the Wilson fits rather than lKin3 jl966l) . 
Points that scatter to the upper-right quadrant of this plot thus rep- 
resent clusters for which the fwilsoni model is the best match of the 
three. Once again, 90% (111/124) of GCs in NGC 5128 have both 
Ak62 > —0.2 and ApL > —0.2, such that IWilsod models fit them 
at least as well as either of the other two models. This is again in 
keeping with the situation for both old and young massive clusters 
in the Mil ky Way, the Magellanic Clouds, and the Fornax dwarf 
spheroidal: iMcLaugh lin & van derMarellJioO^ show that the cur- 
vature inherent in|wilson model envelopes is never any worse, and 
often much better, than an unlimited power law as a description of 
any observed extended halo structure. 

4.5 Model dependence of cluster properties 

Having argued that isotropic lwilsoni models fit the majority of GCs 
as well as or better than any other model that we have tried, we 
next ask whether estimates of physical cluster parameters can be 
significantly biased by fitting the "wrong" type of model. 



Table 10. NGC 5128 GCs with poorly constrained core parameters " 



Name 




Name 




Name 




(1) 


(2) 


(1) 


(2) 


(1) 


(2) 


AAT115339 


+4.70 


C032 


-1.37 


G284 


-5.11 


AAT118198 


-2.60 


C116 


+4.54 


K131 


-3.41 


C003 


-0.94 


C129 


+4.82 


PFF079 


-3.53 


C03() 


-1.03 


F2GC69 


-7.43 


R223 


-4.99 



" (5/iv,o = /Jv,o(Wilson) -/j,,o(King 1966) 



Figure [To] compares values from lWilsod and lKingI ( 1 19661) fits 
for the central surface brightness /Ji/o, projected core or half-power 
radius log R^, central velocity dispersion log CTp o, effective radius 
log Rt,, total luminosity log Ly, and global binding energy log E\,. 
The values of all parameters are taken from Tables |7] and [8] above, 
and their differences between the two models are plotted against 
the relative quality-of-fit parameter A defined in equation l ll4t . We 
again do this for only 124 clusters; C168, F1GC34, F2GC14, C145, 
C152, C156, and C177 are excluded as before. 

The open circles in all panels of Figure [10] refer to clusters 
with differences of more than 0.6 mag arcsec"- between the fitted 
yUvo in the two models. There are 12 of these, listed for reference in 
Table [To] Arrows are attached to the points for them in the figure 
when the parameter differences fall outside the plotted range in any 
panel. The innermost structures (/Jv.o, ^c, and CTp o) of these clusters 
are not well constrained in a model-independent sense, even if more 
global properties like /?;,, Ly, and E], are generally better behaved. 

These 12 objects aside, the two models return encourag- 
ingly consistent cluster parameters. The mean and rms scatter 
of the difference ^/j^.o = [/iv,o(Wilson) - //v/o(King 1966)] is 
0.10 ± 0.15 mag arcsec"' (compared with an rms errorbar of 
± 0.3 mag arcsec"^), while the average <(5(log R^Yl = 0.05 ± 
0.05 dex (rms) versus a typical errorbar of ±0. 12. Thu s, in or- 
der to accomodate the more distended haloes of IWilsonl models, 
th e data force these fits to (slightly) larger and fainter cores than 
m [Kind ( 1 19661) models. These effects combine to make the aver- 
age ((5(log CTpo)) ~ 0.5 <5(log Re) - 0.2(5/ii/o) a com pletely negligi- 
ble 0.004 dex. Similarly, the average Wilson-minus- lKing offset in 
global binding energy (which in either model is dominated by the 
core) is only ((5(log Et)) - 0.02 + 0.02 dex (rms), compared to an 
rms errorbar of ± 0. 10 dex. 

On the other hand, the projected half-light radius and total 
luminosity differences are slightly larger, and more systematic: 
<(5(log ^,,)> = 0.07 ± 0.08 dex and <5(log Ly)) = 0.05 ± 0.04 dex, 
compared to average errorbars of ± 0.04 dex in both cases. The lu- 
minosity offse t is the same as the ^ 0.15-mag brighter integrated 
magnitudes of IWilsoj models, which was seen in Figure [7] above. 
This and the systemati cally larg er Ri, both result fr om the intrinsi- 
cally stronger haloes of lWilson versus [Kind d 19661) models. 

In most of these examples, the differences between the lWilsor] 
and King ( 1966) parameter estimates tend to be slightly larger for 
more negative values of the relative statisti c A, i.e., for clusters 
that are significantly better fit by the IWilsonl sphere. Clearly, the 
parameter values from this fit are then to be considered the more 
reliable of the two. But even so, we conclude from this discussion 
that analyses o f structural c orrelations w ill not b e significantly af- 
fected by using lKin3 ( 1 19661) models when lwilsoni spheres are better 
fits — or vice versa. The net effect of such an "error" on any given 
correlation will generally be at the level of hundredths of a dex, 
which is within most observational uncertainties. 
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Figure 10. Left panels: Compaiison of isotropic IWilsonI and I Kind ) 19661) model values for the intiinsic central suif ace brightnesses, projected core radii, 
central line-of-sight velocity dispersions, projected half-light (effective) radii, total luminosities, and global binding energies for 124 reliable GC candidates in 
NGC 5128. Open circles represent 12 clusters with |pv/o(Wilson) -//v,o(King)| > 0.6 mag arcsec"^, which have particulai'ly poorly constrained core stmctures. 
These objects ai'e listed in Table llOl Filled points represent the remaining 112 clusters whose core parameters are better constrained. 
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Figure 11. Similar to Figure [To] but comparing cluster properties derived from ISersl3 versus IWilsoiJ model fits. Open stars represent 18 clusters with 
|//v/o(Sersic) - yUv.o(Wilson)] < -4.0 mag arcsec"^. 
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Figure 1111 compares estimates of the same cluster parameters 
from lSersid models versus IwilsonI models. The open stars in these 
graphs refer to 18 clusters with particularly large |jUv/o(Sersic) - 
^v,o(Wilson)] < -4 mag arcsec~ '. Even apart from these, it is 
immediately obvious that ISersid models typically imply much 
brighter and sm aller cor es, and rather low central velocity disper- 
sions, relative to lWilsor] (or, therefore, lKinslll966h models. This is 
mainly a combination of model artifact and the PSF-blurring of the 
cores of our clusters. 

With an assumed I{R) ~ exp(-i?'^"). ISersicl models are pushed 
towards high n to fit the extended haloes of the observed GCs: 
more than half (69/124) of the clusters represented in Figure [TT] 
have « > 2, which implies very strongly peaked central density 
profiles even in projection (e.g., see Figure |5j also, recall that the 
unprojected density diverges and the velocity dispersion vanishes 
at r = for any n > 1). These are probably allowed by the data 
here only because the PSF flattens the models; realistically, most 
glob ular clu sters do not have such steep and non-isothermal cores0 

ISersid models are therefore not a good choice to characterise 
glo bular c lusters whose cores are not fully resolved. Nevertheless, 
the lSersi3 values of the more globally oriented structural parame- 
ters in Figure [TT] (/?;,, Ly, and £;,) are in better overall agreement 
with those returned bv lWilsonl and lKind ( 1 19661) model fits. These 
quantities are evidently rather robust, in general, against the de- 
tails ofttiemodel adopted to measure them. All in all, we conclude 
that lWilsonl fits can safely be used to investigate any GC parame- 
ter correlations, and it is indeed the model th at we emp loy for our 
discussion of the cluster fundamental plane in lPaper IVl 

4.6 Clusters observed in lPaper j 

To build a cluster database for NGC 5128 that is as complete 
as practic al, we would like to fold the sample of 27 GCs ob- 
served in IPaper j jHarris etal.ll2002h into o ur curr ent analysis^ 
As we mentioned at the beginning of ^ iKind ( [T966i) mod- 
els were fitted to STIS and WFPC2 surface-brightness profiles 
for t hese clusters. However, a different modeling code was used 
( [Hol land. Cote, & Hesser 1999) to obtain a smaller set of cluster 
parameters than we have derived for the ACS sample here. We do 
not re-fit the original data, especially as 9 of the clusters there have 
in fact been re-observed and re-modeled as part of this paper. In- 
stead, we calculate the full suite of physical properties in Tables[7] 
[8] and [9] within the context of King (1966) models only, starting 
from the few parameters given in Table 2 of jPaper j . 

It is necessary first to correct all of the intrinsic central sur- 
face brightnesses pre viously published for these GCs . The earlier 
model-fitting code of [Holland. Cote. & Hessei^ ( [l999h did not ex- 
plicitly compute yUv/o for t he best King ( 1966) model, so it was es- 
timated for every object in [Paper j as the directly measured bright- 

^ We have, in fact, re-fit th e 153 non-c ore-collapsed clusters in the Lo- 
cal Group sample of McLa ughlin & va n der Maiel ( 2005) — the cores of 
which are all completely resolved — with Sersic models. Only 22% of these 
(34/153) have ii > 2; and regardless o f the va lue of n, only about 5% of all 
clusters are sig nificantly better fit by a iSersicj model than by the isothermal 
cores of lKindfT96ffl or , Wilson models. 

An entirely d ifferent samp le of 21 clus ters was also observed wi t h HST , 
and fitted with [Kind 119661) models, by [Holland. Cote, & HesseJ h999[) . 
However, their imaging program used the lower-resolution WF chips of the 
WFPC2 camera and focused on clusters at smaller galactocentric radius, 
against higher sky levels, than almost all of ours. We have not attempted to 
incorporate any results from these lower-quality data into our work. 



ness of the central 1 pixel of the cluster images. But the STIS pixel 
scale is 0'/0508 px"', so that 1 px = 0.94 pc. Unfortunately, this is 
comparable both to the half-width at half-maximum of the stellar 
PSF and to the typical GC core radius of Rc ~ 1 pc. The intensity 
of the innermost pixel on a cluster image in this case is therefore 
an average of the intrinsic surface-brightness profile over the en- 
tire core and beyond, which is much fainter than the brightness at 

exactly R = 0. 

Here we use numerical calculations of [Kind ( Il966h models 
to derive a self-consistent fiv,o from the published Wp, angular 
ro, and total apparent magnitudes Vtoi of the clusters in [Paper it 
y"v,o = V'tot - 2.5 log [lo/iJlrl)] for ro in arcsec, where Iq and X are 
dimensionless functions of Wq (o r the concent ration parameter c) 
described in, e.g.. Appendix B o f [McLaughlii] jloOQ) . Extinction- 
corrected Vtot are o btained fr om the Washington CTi photometry 
given m Table 3 oflPaperll' and our Table [4| here, in exactly the 
same way as described in connection with Figure [7] above. We as- 
sign an errorbar of ± 0.3 to every Wo value, ± 10% to every ro, and 
±0.05 mag to every V,„t, and propagate these to obtain uncertain- 
ties for the re-derived ^vo. Then, using the projecte d galactocentric 
radii given for all the clusters in Table 1 of Paper I and population- 
synthesis mass-to-light ratios calculated for them in Table[5]above, 
we derive all of the parameters that we did for our ACS cluster 
sample in 94.21 

Tables [AT[lA4l in Appendix [Alcontain the results of this. Like 
Tables (eJlD they are available in machine-readable format, either 
online or upon request from the first author. 

As expected, the fiv,o published in Table 2 of [Paper 3 are all 
significantly fainter than our corre cted val ues in Column (11) of 
Table I All even after correcting the [PaperT numbers for extinction 
(Ay = 0.341 mag in all cases), we still have an average difference 
of (yUv,o(new) - ;Ui/o(old)> = -0.72 ± 0.05 mag arcsec"^, due to the 
uncorrected effects of seeing on the old estimates. 

Note that our effective radii Rj, in Column (6) of Table lA2l are 
smaller than the quantities r;, listed in Table 2 of Paperj^his is 
because we record the projected half-light radii, while [Paper j re- 
ported depr ojected value s. The ratio RiJri, ranges between ^ 0.73- 
0.76 for the [Kind d 1 9661) models that appl y to real clusters. 

The 9 clusters in common between ' Paper ll and the current 
ACS sample (C007, C025, C029, C032, C037, CI04, C105, G221, 
and G293) are included in Tables rAT[4A4l We have compared the 
basic parameters c, Vtot, rQ, and Ri, g iven there to the results, in Ta- 
bles|6]and|7| of our new [Kind d 1 9661) model fits to the ACS profiles 
of these objects. The old and new concentrations and the integrated 
cluster magnitudes agree quite well, but the radii ro and Rj, do not: 
both scales are on average smaller, by a factor of about 0.7 (0.15 
dex) in the mean, in our new analysis. As a result of this, the intrin- 
sic central surface brightnesses for these 9 objects in Table [6| are 
typically brighter (by 0.8 V mag arcsec"' on average ) than even 
the revised estimates in Table lAll based on the [Paper U rn values. 

We have not been able to identify the cause of this discrepancy. 
Since the very well resolv ed Ri, di ffer by as much as the core scales 
ro between this paper and [Paperl ! for these 9 GCs, our PSF model- 
ing and convolution are not responsible. Also, direct inspection of 
the ACS data shows clearly that the clusters are not as intri nsically 
faint at thei r centres as t he fiv,o derived from the large ro in [Paper j 
imply. Our [Kind ( [19661) fits in ^ above should be take n as the 
definitive ones for the clusters observed both here and in Paper j 
For the other 18 GCs in the earlier sample, we caution that any 
structural or dynamical properties derived through ro or Ri, could 
potentially be biased by overestimation of these radii in [Paper 3 . 
This applies to most of the parameters tabulated in Appendix [A] 
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Table 11. Predicted aperture velocity dispersions from 147 profiles of 131 GCs in NGC 5128 



(1) 



[Mo L-^y\ 



(2) 



(4) 



log Rh 
[pc] 
(5) 



log Ri, 
[arcsec] 
(6) 



log (ri,p(«,,/8) 
(7) 



log i^ap(«/,/4) 
[kms-'] 
(8) 



log 0-!ip(«,|) 

(9) 



log o-ap(4R,,) 
(10) 



log crap (8 K,,) 
[km s-'] 
(H) 



AATl 1 1563 WFC/F606 



AATl 13992 WFC/F606 2.914;];*^5 



AATl 15339 WFC/F606 2.015^ 



K66 

W 

S 

K66 

W 

S 

K66 

W 

S 



0.468 
0.486 
0.478 
0.560 
0.749 
0.528 
0.376 
0.397 
0.384 



_0798+0.005 
-0.780!8f 

-0 705+''*'' 

-0.516„|„, 

n -7TO+0.016 
-0.738„„j5 

-0.890+''™' 

-0.869+||-i|?5 



0.835+j!- 
0.838+ij- 
0.812+!!' 
0.893+°- 
0.893+jj- 
0.863+J 
1.034+9- 
0.970+°- 
0.953+°- 



.027 

e9 



n 833+0-027 
U.8J3 

U.533_(,„3|, 



879+0.032 
"-""-0.035 

1 022+0-"27 

Urn 



0.967 
0.967 



A o 1 I +0.027 

am 

-0.030 



0.809+ 



A s I ^+o:027 
-9.Q3Q 

0.856 ujj, 
8 71 +0.032 
''■'*"-0.036 
0.977+0.027 

Q 94Q+8:83? 
"-0.029 



0.950+ 



0.762;"- 
0.762;°- 
0.764 " 
0.820_ 
0.807^^^ 
0.822;°- 
0.914+I!- 
0.893 
0.898 



:815 



.038 
.032 
.036 
.026 

mi 



Q752+0.027 

Q 749+8:859 
'^^-0.030 
Q754+0.027 

Q g I ^+8:8^9 

789+'^-^^^ 
"-"'^-0.038 
n 8 1 fi+0.032 
"-'"^-0.036 
Q9Q7+0.026 

879+'^-^^^ 
-0.029 
n 887+0-0-^' 



A machine-readable version of the full Tablelslis available online Chttp://vi'ww.astro.keele.ac.uk/~dem/clusters.html) or upon request from the first author. Only a short extract from it is 
shown here, for guidance regaiding its form and content. 



and, e.g., to dy namical masses estimated from velocity- dispersion 
measurements jMartini & Hdl2004l ; lReikuba et al.ll2007h . 



5 PREDICTED VELOCITY DISPERSIONS 

Apart from studies of parameter correlations, we expect that one 
eventual use for the material presented in ^ will be to facili- 
tate the comparison of dynamical GC mass-to-light ratios (derived 
from internal velocity-dispersion measurements) to the population- 
synthesis model values that we calculated in 93.11 Such compar- 
isons could potentially be used to constrain some of the stellar- 
population aspects of these models; perhaps to identify clusters 
much younger than the 13 Gyr we generally assume; and also to 
examine the extent to which long-term dynamical cluster evolution 
(e.g., the evaporation of low-mass stars driven by two-body relax- 
ation) might affect M/L. 

We have therefore calculated aperture velocity-dispersion pro- 
files for all 131 objects in our ACS sample using the best-fit struc- 
tural models from each of the families we have explored in detail, 
and the 13-Gyr T^°'' listed in Table|5] Given any of our models with 
a fitted yq and Wq (or c) or Sersic ;i, solving Poisson's and Jeans' 
equations and projecting along the line of sight yields a dimension- 
less CT'p = iTp(i?)/a"p,o as a function of projected clustercentric ra- 
dius R. The predicted velocity dispersion within a circular aperture 
of any radius /?ap then follows from 



1 20071) have done so for 27 clusters including some of the 



p.o 



J-^Rap 




all{R)RdR 



I(R)RdR 



(15) 



Using CTp from Table [8] we have obtained cTap for every 
structural model for every cluster, within each of five apertures: 
^ap = RiJi, Ri,/4, Rh, ARh, and 8i?;,. The results are listed in Ta- 
ble | 1 1 1 Including cTp o itself gives cr^^ at = for i Kind jl966h and 



1 

IWilson models, or ati^ap = Rc f or l S ersid models . The dispersion in- 
side any other aperture can then be obtained from interpolation on 
our numbers. 

By far the majority of a cluster is contained within 8 Ri, in any 
realistic model. Thus, the velocity dispersi on in any l arger aperture 
diflfers from crap(8i;ft) by <0.01 d ex for a ll lKhi^ fl966) and .WilsonI 
models, and by < 0.02 dex for all lSersij models with n < 6 (which 
includes about 95% of the GCs here). In any case, a typical GC 
half-light radius of 3-4 pc translates to a'16-a'22 in NGC 5128. 
A slit of width 1", for example, therefore usually corresponds to 
if ap I Rh of Just a few . 

iMartini & Hoi ( |2004|) have measured the veloc ity dispersions 
for 14 relatively bright GCs in NGC 5128, and iReikuba et al.1 



iMartin i & Ho' objects. The combined sample from these two stud- 
ies includes 8 globulars modeled in this paper that d o not ap pear 
in Table [TO] above, and 9 objects obser ved only in IPaper j . for 
which we have re-computed iKinj ( 1 19661) parameters in Appendix 
lAl Thus, we no w compare the observed (Tap of these 17 GCs to 
our lKin3 ( fl96^ model predictions for them. (We have calculated, 
though not tabulated, model velocity-dispersion profiles for the 
IPaper j clusters in exactly the same way as for the current ACS 
sample.) 

For clus t ers ob served by both IMartini & Hoi ( |2004|) and 
iReikuba et al] ( |2007|) . we take the root-mean-square of their in- 
dependent velocity-dispersion measurements to est imate a single 
observed cTap. Then, since both IMartini & Hcl and IReikuba et al.l 
used slits of width 1" for their observations, we adopt an effec- 
tive aperture radius of R.^^ = 0'.'71 — half the diagonal length of 
a 1" square — to specify a predicted (Tap for every cluster. With 
D = 3.8 Mpc for NGC 5128, this corresponds to /?ao ^ 13 .0 pc, 
which is significantly larger than the projected ( King [196^ half- 
light radii of all 17 GCs being considered here: 1.7 < R^p/Ri, < 5.7, 
with a median of 2.8. As a result, even though the choice of R.^^ is 
not particularly rigorous, changing it by a factor of two affects our 
comparison of dynamical and population-synthesis mass-to-light 
ratios at only the ^5% level. 

The upper panel of Figure [12] shows the observed versus 
model iTap for the 17 objects in question. The lower panel then 
shows the square of the ratio of observed to predicted dispersions, 
'-'"obs/''"pred' ^hich is equal to the ratio of dynamical to population- 
synthesis Yv- The filled circles in both panels of the figure refer to 
clusters with ACS data anal ysed in this paper; the open circles, to 
clusters only in the IPaperl sample. The apparent mean off^set be- 
tween these two sets of points could be the result o f overes timating 
the half-light radii of the GCs measured only in IPaper ll (as dis- 
cussed at the end of 34.6l l, since a spuriously large Ri, implies a low 
'-'"pred '^^'^ Lv I Rh- In any event, the median ratio of tTobs/Cpred for 
all 17 clusters is 0.90, and the rms scatter about this value is ± 0.20. 
The median and rms scatter of Y^^"/Ty°'' is therefore 0.81 ± 0.40. 

The main point to be taken from this is that our population- 
synthesis estimates are relatively good guides, on average, to the 
global dynamical mass-to-light ratios of globular clusters in NGC 
5128. This is also true in general of GCs in the Local Group that 
have been modeled in the s ame way (McLaughlin & van der Mar^ 
2005; Bar mbvetal.ll2007i). In our fun damental-plane analysis of 
iPaper IV (and also lBarmbv et al.ll2007t) . we therefore use the clus- 
ter properties derived from T^"*" in 94.21 to explore accurately a 
wide range of structural and dynamical correlations for much larger 
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Figure 12. Top/ Velo city dispersions ob s erved for 17 GCs in NGC 5128 by 
iMartini & Hoi j2004) and lReikuba et alj <2007l) . a gainst dispersions within 
apertures of radius ff.'ll = 13 pc as predicted by our iKind il966h model 
fits and population-synthesis mass-to-light ratios. Bottom: The implied ra- 
tio of dynamical and population-synthesis mass-to-light ratios, as a func- 
tion of model T^"''. In both panels, filled circles denote 8 clusters in the 
ACS sample of this paper: C006, C007, C012, C018, C025, C()29, C036, 
and C037. Open circles correspond to 9 GCs observed only in Paper I, 
with revised King ( 1966) structural parameters given in Appendix IaI here : 
C002, coil, C017, C021, C022, C023, C031, C041, and C044. The bold, 
dash-dot lines indicate the median ratios o"ap(obs)/(Tiip(pred) = 0.90 and 
T^/7t7 = 0.81. 



samples of globulars than would be possible if we were restricted 
to objects with direct velocity-dispersion measurements. 

A more specific imp lication is that, as br ight and as massive as 
the GCs studied by Mart ini & Hoi ( |2004|) and lReikuba et al.l ilOOH) 
are (the 17 in Figure [12] have 6.3 x 10^ Lq < Ly < 2.6 x 10** Lq 
and an average (Tv> ~ 2 Mq L^' in our analysis), there is noth- 
ing in their mass-to-ligh t ratios alone to s uggest that th ey are 
otherwise unusual. Both iMartini & Hoi and iRejkuba et alj quote 
an average dynamical (T^*") 3 Mq Lg' for their cluster sam- 
ples, which is a factor of two higher than the average for a rep- 
resentative sample of Milky Way GCs (e.g. McLaughlin 2000; 
iMcLaughlin & van der MarelllioOSi) . But nearly half of the NGC 
5128 clusters plotted in Figure [T2] are expected to have Tv in this 
range. This is simply because they are redder and, by inference, 
more metal-rich than the majority of Galactic globulars; see Fig- 



ure [3] above, and note also that T!^"'' along the x-axis of the lower 
panel in Figure [T2| is a rough placeholder for cluster metallicity or 
Washington (C - T\)q colour. The other half of the GCs in Fig- 
ure [T2| are relatively bluer and more metal-poor, and they have an 
expected T^"'' a 2 Mq Lg' that is closer to the Milky Way GCs. 
While 4 of these appear to have larger true ,tJ^", 5 others do not. 
Moreover , the clu sters with the largest Y^^"/Ty°'' appear only in the 



sample of Paper j. which again may have overestimated Ri, and thus 

Aside from any possible issues with Ri, measurements, many 
of the dynamical mass-to-light ratios that we have obtained for the 
17 GCs in Figure 1 121 a re smaller than those computed by either 
IMartini & HQl hooiT or lReikuba et alj ( |2007|) for the same objects. 
This is usually because both of these other groups estimate "virial" 
cluster masses using the formula M^-„ = 10 <rl^Ri,/G, where cr„ is 
meant to be the velocity dispersion measured through an infinite 
aperture and Ri, is the projected half-light radius. However, the co- 
effi cient in this expression is model-dependent, ranging from 8-10 
for iKing' ('1966') models with concentrations c < 2.5. In addition, 
both|Martini & Ho and Re ikuba et alj assume that the aperture cor- 
rection to relate their measured (T abs to the formal (Too is always 
negligible. In fact, our iKingj l l 19661) model calculations for the clus- 
ters in Figure [T2| yield cr^/(T^p(0'.'71) 0.75-0.96, so making the 
correction (which we effectively do) can sometimes lower Mvi, and 
T^^" noticeably. 

Finally, we note that the ratio Y'^^/Y'y^ does not correlate 
with the total luminosity of the clusters in Figure [12] We conclude 
again that these objects a re not special or unusu al in terms of their 
mass-to-light ratios only. Martini & Hoi ( |2004|) and Reiku ba et al.l 
l i2007) go on to consider how the massive GCs in NGC 5128 relate 
to Milky Way globulars, dwarf-galaxy nuclei, and ultra-compact 
dwarf galaxies, in terms of structural correlations and the funda- 
mental plane. We will discuss this further in .Paper IV, . 



6 SUMMARY 

We have used new HST imaging with the ACS/WFC camera to 
construct internal surface-brightness profiles and derive stmctural 
parameters for a sample of 131 globular cluster (GC) candidates in 
NGC 5128, the nearest giant elliptical galaxy. We combined these 
with a s ample of 27 clust er s previo usly measured with STIS or 
WFPC2 jHarris et al.l2002l , = |Paper t . for which we have corrected 
and extended the structural parameters given in Paper j 

After identifying 4 ACS-observed objects that may not be 
not bona fide GCs and another 3 that have especially ragged or 
contaminated surface-brightness profiles, and accounting for 9 re- 
observations from the Paper I sample, we have now observed and 
modeled a total of 142 distinct, reliable GC candidates in NGC 
5128. The catalogue of para meters that w e have produced for 
them rivals the compilation o f lHarrisI ( 1 19961) for Milky Way glob- 
ulars in size and is based on more homogeneous data and analy- 
sis. The final product is very simila r in most respects to thos e of 
iMcLaughlin & van der MardI l l2005h and lBarmbv et al.l ( |2007|) for 
a total of 196 old GCs and 50 massive young star clusters in the 
Milky Way, the Magellanic Clouds, the Fornax dwarf spheroidal, 
andM31. 

We fitted the NGC 5128 clusters' surface-brightness profiles 
with three different structural models: the standard King ( 196g) 
modified isothermal sphere, an isotropic version of lwilsoa 's tl97a) 
alternate modificat ion of an isoth ermal sphere, and the surface- 
density profile of ISersij ( Il968l) . We estimated K-band mass-to- 
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light ratios individually for every cluster, by applying population- 
synthesis models given estimates of cluster metallicities from their 
(C-T\) colours in the Washington filter system and assuming a uni- 
formly old age of 13 Gyr. We combined these with our structural 
modeling to derive a wide range of physical properties including 
central surface brightness and central potential, concentration in- 
dices, core and half-light (effective) radii, total luminosity, and pre- 
dicted dynamical parameters including total mass, binding energy, 
central mass density, central escape velocity, relaxation timescales 
and phase-space densities, and velocity dispersions in series of cir- 
cular apertures. 

Extensive consis tency checks and i nter compar isons of the 
three different models ( King 1966. Wilsonl and lSersi3) have shown 
that in most cases the Wilson model fits at least as well as the oth- 
ers and in many cases is superior to them, taking better account of 
the detailed shape of the distended, low-intensity outer envelopes 
that many clusters have. In this respect the NGC 5128 globulars 
are like those in the Mil ky Way, the Magellanic Clouds, and the 
Fornax dwarf spheroidal jMcLaughlin & van der Marelll20 05l), al- 
though the situation in M31 is not as clear jBarmbv et al.| [2007). 

We compared aperture velocity dispersions predicted by our 
analysis for 17 GCs in NGC 5128, to actual measurements by 
iMartini & Hoi ( l2004h and lReikuba et al.l ( l2007t) . The expectations 
based on population-synthesis models are in reasonable overall 
agree ment with the observations. 

In lPaper IVl jMcLaughlin et alj2 Q07'). we use the data derived 
here to discuss a wide range of structural correlations for globular 
clusters in NGC 5 128, including comparisons with Milky Way GCs 
and w ith other types of massive, cluster-like objects. ISarmbv et al.l 
| |2007|) also use our results, in a discussion of the combined funda- 
mental plane of old GCs in six galaxies. 
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APPENDIX A: PARAMETERS OF CLUSTERS 
OBSERVED IN PAPER I 

Here we pr esent the tables of fitted and derived properties for 
iKind ( Il966h model descriptions of the 27 globular clus ters ob- 
served by STIS o r WFPC2 and previously modeled in IPaper ll 
( [Harris etalj|2002l) . The construction of these tables is detailed 
in 94.61 Their contents and format are essentially the same as 
those of Tables |6H9] for our main sample of clusters observed 
with ACS (see the descriptions in 34.11 and 34. 2K and they 
are also available in machine-readable format, either online (at 
http://www.astro.kee le.ac.uk/~dern/clusters.html) or upon request 
from the first author. 

Some of the columns in Tables IAmA4l are present only to 
keep all of our tabulations of the fits to the combined ACS-t jPaper j 
cluster sample in a single, well-defined format. For example. Col- 
umn (4) of Table lAll contains the "colour" V - V, which is of 
course always identically zero but is the analogue of the nontriv- 
ial {V - F606)o in Table |6l(M above). 



© 2007 RAS, MNRAS 000.[Tlt26l 



Globular clusters in NGC 5128 



Table Al. Basic parameters of lKind fT96^ model fits to 27 GCs from lHarris et alj j2002l = |Paper ll) 



Name 


Detector 


Ay 


V -V 


/Vpls 


Model 




^bkg 






[mag] 


[mag] 








[Lo pc-2] 


CD 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


C002 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C007 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


COl 1 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C017 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C021 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C022 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C023 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C025 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C029 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C031 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C032 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C037 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C040 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C041 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C044 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


ClOO 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


ClOl 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C102 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C103 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C104 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C105 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


C106 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


G019 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


G22! 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


G277 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


G293 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 


G302 


HAR02/V 


0.341 


0.000 ± 0.000 


5 


K66 


1.00 


0.00 ± 0.00 



Wo 



(9) 



(10) 



[mag arcsec^' 
(II) 



log ro 
[arcsec] 
(12) 



log ro 
[pc] 
(13) 



8.50!g;» 
8.00;0;» 
o 90+0.30 

6 60+° ™ 

-0.30 

o in+0.30 

*'-'"-0.30 

7 30+° ™ 
'■^"-0.30 
7 50+° ™ 
'■■'"-0.30 
g 20+0.30 

g [0+8:^8 

7 70+'^-^'^ 
'■'"-0.30 
n+0.30 



i.SO^j 
5. 10+! 



-0.30 
75Q+0.30 

^■"'-0.30 
70o+°-30 

7 7<;+8:][] 
-0.30 
55Q+0.30 

9 20+'^-^'^ 
-0.30 
,n+0.30 



7.60^! 
7.20+! 



ft on+o.3() 

74Q+0.30 

7 10+''-^'' 
■ -0.30 

6 90+°-™ 
"■'"-0.30 

7 sn+°-30 
'■''"-0.30 
7 I <;+0.30 
'■^■^-0.30 



1 89+0.09 
'■=^-0.01) 
1 89+°-»'' 



, ^9+0.09 
'■'52-0.01) 
I ftS+0.09 
'■6»_0.09 



1 74+0.09 
'■'^-0.09 

2 07+0.08 

I Sfi+^-8^ 
'■''''-0.09 
^ 68+° "^ 

L86+E 



1 65+°"'' 

'■-°-0.09 

1 50+°"'' 
'■^"-0.09 
I 77+0.09 
-0.09 
1 ^7+0.09 
'■- '-0.09 



1478 



15.31 
14.91 
14.57 
13.26; 
15.19; 
15.10" 
14.51 
14.21 
13.90; 
16.80; 
14.64 
15.15; 
16.92" 
18.03 
19.07 
15.35 
16.09; 
21.24" 
16.27 
16.70 
15.70; 
15.59; 
14.70" 



,+0.08 
'-0.10 
I +0.10 
-0.11 
1+0.09 
-0.10 
+0.14 


1 ^/;-7 +0.041 
_] 119+0.041 


-0 101+"°'" 
"■'"'-0.046 
0145+0.041 


1 rcc+0"04I 
'■"■'-0.046 
-0.914+°"^' 


n 1 1 n+8'83l 
"■l "'-0.046 
0352+0.041 


:8:J§ 

-0.11 
f+0.12 
-0.13 
-+0.12 
-0.13 
1+0.09 


1 1 ci-r+O (Ml 

-1.229+°"'' 

I T^n+O^fwY 

-'■328S 
-1.276+°"" 


0.078;»;g« 
0.036+"°'" 
-0.063!^« 
-0.010+"°'" 




1 ) 94+8:8^1" 


O.U/Z 


-0.12 
+0 07 


-1 310+''-''^^ 
'■^'"-0.046 

1 CT>+O041 

— 1 .523 


-0 044+''-''^^ 
"■"^-0.046 
t\ 1CO+0 041 


,m 

'-0.]1 

)+9-!l 


-1 509+''-^^ 
'■^"^-0.046 

-0.967+J"! 


-0 243+°^°'' 
"■^*^-0.046 

0.299+J9^! 


,:8:i^ 

-0.11 
;+0.13 

,:8:lf 

■-0.12 
(+0.14 


1 977+8:8ti 

'■"'-0.046 
174+0.041 

-1 066+''-^^ 
' ■"" -0.046 
-0 959+°"4' 


-0 111+"MI 
"■'"-0.046 
0.091+"°'" 

200+^ '^^^ 
"■^""-0.046 
0307+0.041 


,:8:J^ 

-0.09 
;+0.]l 

,:8:I5 

-0.14 
1+0.12 

1:8:11 

-0.15 
1+0.12 


1 999+8:811' 
'■^^■'-0.046 
-1.27S+°"*' 


043+^-'^^^ 
"■"^^-0.046 

-0.010+"°'" 


-1 252+''-^^ 
' ■^■'■'-0.046 

-0.719+°"" 


014+8:815 

"■"'^-0.046 
0545+0.041 


_14jj+8:8Jt 
'■^■'°-0.046 
-0.99S+°^*" 


-0 191+8:815 

"■"'-0.046 
0.270+"°'" 


:8:l] 

-0.14 
1+0. 13 


1 999+8:8tt 

'■■'^''-0.046 
174+0.041 


043+8:'^^^ 
"■"^^-0.046 

0.091+"°'" 


,:8:I5 

-0.12 
1+0.13 
-0.14 


_ 1 444+8:81? 
'■** -0.046 
1 9fis+°041 

-l-268_o„4j 


A 178+8:815 

""■'■"'-0.046 

-0 002+"°'" 
"■""^-0.046 



A machine-readable version of Table fAl l is available online (http://www.astro.keele.ac.uk/-dem/clusters.html) or upon request from the first author. 

Table A2. Derived structural and photometric parameters from lKiri3 il96d) model fits to 27 GCs from lHarris et alj i2002l = |Paper ll) 



Name Detector Model log r^\^^ log Rc log Rj, 

[pc] [pc] [pc] 

(1) (2) (3) (4) (5) (6) 



log («,,/«,) 



log la 
[Ley pc"- 
(8) 



log 70 
[i-G.VPc""'] 
(9) 



log Lv 

[io.l'] 
(10) 



[mag 
(11) 



log Ih _ 
[Le v pc"" 
(12) 



C002 
C007 

con 

C017 
C021 
C022 
C023 
C025 
C029 
C031 
C032 
C037 
C040 
C041 
C044 
ClOO 
ClOl 
C102 
C103 
C104 
C105 
C106 
G019 
022 1 
G277 
G293 
G302 



HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 
HAR02/V 



K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 
K66 



1.88+ 



Q+0.04 
^-0.05 

1 98+"°-' 

■ -Q.Q5 

2 00+"^°-' 

■ -Q.Q5 
I '7^;+0.04 
'■■^^-0.04 
1 94+0.04 

1.65 



).05 
).05 

:8:8? 



,:8:85 



1.98«;- 

1 75+"°-l 
■ -0.05 
1 52+0.04 

1 96+'^:''^ 
'■'°-g.05 

1 69+"°'' 
'■"''-0.04 
9 9 I +0.03 
^■^'-0.04 
170+0.05 

1 m+^^i 

'■""-0.05 

2 13+0.05 



1 92+0.09 
'■'^-0.05 

1 60+"°5 
■ -Q.Q5 
1 59+0.04 

1 59+0.05 

1 57+8:8^ 



I +0.040 
' -0.044 
1 33+0.039 

0099+8:8it 

n 99^^+8:8^^^ 
0.325„|„, 

066+°"3'' 
"■" -0.044 
Q+0.038 
"-0.043 

-0 079+°"''' 
"■"' -0.043 



-0.111 + 



0.018+ 



059+°^"3'' 
"■"■^'-0.044 

-0.060+°"''' 

-0 266+'^:^'' 
"■^°°-0.044 

-0.256;°"-''' 
n 9o9+8:8t^ 
"■^"'^-0.043 

-0 124+°^03'' 
"■'^^-0.044 
070+°0-« 

"■is^ii 

0.280+°:S^ 

037+°0*0 
"■"^'-0.045 

-0 026+°^"'' 
"■" -0.044 

-0 006+°"3* 
"■"" -0.043 
n C97+0.038 
"■■^^'-0.043 

-0.214;';°^"3* 

252+'^:^^ 
"■^■^^-0.043 



_0 i93+o:039 

-n 091+8:8^ 



779+0.056 
"■""-0.057 
865+0 °''^ 

892+8'^^^ 
742+8:854 
"■'^^-0.007 
0.828!0;°5' 

"■se^ii 

517+0.035 

771+0.053 
0.771 „g^| 



0.821" 
0.588 



0.049 
+0.041 



0.878+ 



724+"°^^ 
"■'^'-0.060 

507+0 °5' 
0.50/ „g^^ 

-0.030 
539+0.051 
"■"^'-0.049 
5 57+0.020 

846+'^:''''^ 
"■"^-0.038 

680+"°"'^ 

1 1=9+8:835 

■ -0.058 
595+0.038 
"■^"^-0.033 

521+0.026 
"■■^■"-0.021 

1 054+0.026 
'■"■'^-0.021 
0.236;"°"' 

824+8:8M 
"■"'''-0.027 
0.530+"9?3 



4g2+o:o44 
n 495+0.024 



890+"'"° 
"■"^"-0.097 
O 739+0.093 

0793+0.1)97 

n /1 1/; +8:8^1 
"■'t "'-0.045 
762+0°''5 
"■' -0.088 
549+0.072 

596+^^^ 
"■■^'"-0.069 
793+0.097 
"■'"^-0.091 

762+"°'^ 
"■' -0.088 
547+0.085 

990+^^^ 
"■'""-0.100 
7 62+0.095 

n =Q/;+8:8^i 
''■5"'-0.069 

762+"°''5 
"■' -0.088 
04g7+0.062 

56i;8:8iS 

"■""' -0.077 

401+0.049 
"■*"' -0.043 

1 199+0.094 
'■' -0.098 
521+0.082 
"■"■^ -0.072 
n «97+0.069 
"■■''" -0.059 
O =97+0.069 
"■■^^'-0.059 
450+0.057 

572+8:^^^^ 
"■■^'^-0.066 

0.506+"9S 



574+0.088 
5i7+8:8S 



4.65+ 



4.54+0.04 

4 43+8:8^ 
'■*^-0.06 

4 60+°0* 
'■""-0.04 
473+0.05 

5 95+8:8^ 

■^■^^^-0.05 
4 49+0.04 



+0.04 
-0.04 
+0.05 



4.52 
4.76_ 

4 87+SH 

-0.03 

5 00+°04 

3 84+8:8! 
^■"'-0.05 

4 7 1 +0.04 
'■"-0.04 
450+0.06 

3 79+8:85 
^■'"-0.04 

3 35+0.06 

9 93+8:85 
-0.03 

4 49+°05 
'■'''-0.04 
4 1 9+° 05 
'■'■'-0.05 
2 06+° 05 
''■"°-0.05 
405+0.06 

9 oo+8:85 

^■"''-0.05 
42g+0.06 



4 33+0.06 
'■"^-0.05 
4 68+° 05 
4 43+o.oS 



4 46+""" 
435+0.09 
4 14+8:8§ 

3 8n+^-^^ 
3.SU_oo9 

4 23+0.09 

-0.08 
441+0.10 
'■"-0.09 

5 03+0 °'' 
■^■"^-0.08 

4 20+0 °' 
'■■'"-0.08 

4 16+0°'' 
'■'°-0.08 
4 51+0.09 



3 95+u.uy 
■ -0.08 
452+0.09 



2 59+0.88 
— ^ -0.07 
4 14+0.09 
'■"-0.08 
9 89+0.10 

I 93+0.10 
■ -0.09 



3 95+010 
^■"■^-0.09 
+0.09 



= 89+0.02 
-'■'"-0.02 
539+8.02 

5 i2+8:8i 
5 1 0+8:85 

''■"'-0.02 
5 08+0.02 
"■'"'-0.02 

5 97+0.02 
■^■"-0.02 

6 34+0 °2 
"■'"-0.02 
= 89+0.02 
^■""^-0.02 
5 99+0.02 

■ "-Q.02 
5 91+0.02 

5 8fi+8:85 
= ■'"'-0.02 
5+0.02 

:85 



18.05+ 



4.57 
4.15 



5.85; 
5.64; 

c 8i+8:o2 

'•'"-0.02 
5 go+0.02 

5 45+8:85 
-^-0.02 
4 99+0.02 
'■""-0.02 

4 62+"°2 
'•"^-0.02 

5 63+" °^ 
'■°=-0.02 

5 30+" °2 
■^■^"-0.02 

4 30+" °2 
'■''"-0.02 
475+0.02 

5 60+^-'^' 
■"■""-0.02 
550+0.02 

= 1+8:85 
5-61-0.02 
5 59+0.02 

5.57" 



-0.02 



5+0.05 
"-0.05 
1575+0.05 

17 43+8:85 

-Mi 



17.52!°g^ 
17.80;°;| 

16 89+0^* 

-0.05 

18 18+0.05 
18.18 

17 74+li.B5 
"■"-0.05 
1794+0.05 

18.08+11 
Ig 1 1+0.05 

I a /;9+8:85 
18.62„„5 

18.19;g:» 

18.23;°;| 

19 11+™ 
'"■"-0.05 

20.24;°;| 
21.18^ 

^-^+0.05 

19 48+^-^^ 
9] q7+0.05 
-0.05 

20.85!°! 
18.72!«;| 
18.98!°;» 

I 8 99+0.05 

"'■'--0.05 
5+0.05 



3.51 + 



3.54+' 



3 90+° 00 
"■""-0.01 

3 63+° 08 
"■'"'-0.08 

4 04+° 03 
'■"'-0.04 
4 50+°04 
'■■^"-0.05 
3 48+° 08 
"■'"-0.09 
3 55+0.08 

-0.08 



+o:o5 



18.75+ 
18.82+! 



3.94 



S^OS-oS 
3 74+0.08 
-0.08 
389+0.01 

2 96+'^-'^^ 
■'■" -0.07 
9 89+0.01 
^■""-0.01 

1 50+°'° 
'■■^"-0.08 
354+0.05 

3 46+''-''^ 
"■'"-0.03 
1 40+°02 
'■'"-0.03 



354+o:o2 

9 99+8:8? 

"■'"-0.01 
3 83+0.06 



A machine-readable version of Table fA2| i.s available online (http://www.astro.keele.ac.uk/-dem/elusters.html) or upon request from the first author. 
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28 D. E. McLaughlin et al. 

Table A3. Derived dynamical parameters from lKin^ )l96d) model fits to 27 GCs from lHarris et alj )2002i = |Paper ll) 



(1) 



[Mo Ll\] 
(3) 



log Mioi 
[Mo] 
(5) 



log 
[erg] 
(6) 



log So 

[Afepc"-] 
(7) 



logpo 
[Mope--'] 
(8) 



log E,, 
[Mo pe-- 
(9) 



log o-p 
[km s-'] 
(10) 



VescO 



log 
[kms-'] 
(11) 



log frh 

[yr] 

(12) 



(13) 



C002 


HAR02/V 


2 Q, 3+0.237 
^■""-0.234 


K66 


C007 


HAR02/V 


1 998+0.237 


K66 


COl 1 


HAR02/V 


K66 


C017 


HAR02/V 


1 904+0.242 
'■^'"-0.241 


K66 


C021 
C022 
C023 


HAR02/V 
HAR02/V 
HAR02/V 


1 ncj+O 240 

r (-1-70+0237 

T -7T7+0 '40] 

' -0.393 


K66 
K66 
K66 


C025 


HAR02/V 


T gi 1+0.430 
^■^^'-0.422 


K66 


C029 


HAR02/V 


2 833 

1 1 1 1+0462 
3. 1 

3064-0.447 


K66 


C031 


HAR02/V 


K66 


C032 


HAR02/V 


K66 


C037 
C040 


HAR02/V 
HAR02/V 


1 -1CO+0 270 

2.253 

2 1 38+83tt 
^■""-0.245 


K66 
K66 


C041 


HAR02/V 


2 990+0.444 
^■'^"-0.436 


K66 


C044 
ClOO 
ClOl 


HAR02/V 
HAR02/V 
HAR02/V 


1 915+0.241 

1 898+8:^^ 

2 409+8:533 


K66 
K66 
K66 


C102 


HAR02/V 


2 217+**-2^' 
3:027+|i 


K66 


C103 


HAR02/V 


K66 


C104 


HAR02/V 


1.896;0;^« 


K66 


C105 


HAR02/V 


7 S77+0.412 
'^■''■^'-0.404 


K66 


C106 
G019 


HAR02/V 
HAR02A' 


2 709+"-:^^:' 
1 904+8:5^^ 

'■^"*-0.241 


K66 
K66 


G221 


HAR02A' 


7 303+0.281 
^■^"^-0.274 


K66 


G277 


HAR02A' 


1 994+0.237 
' -''^-0.234 


K66 


G293 


HAR02/V 


1 OO7+0.248 
l-»°^_0.248 


K66 


G302 


HAR02/V 


1 921+0.235 
-0.238 


K66 



6.18 
6.69 



+0.05 
-0.06 
+0.05 



Ik f. I +0.06 
6.61 

54Q+0.05 
527+8:8§ 

^ 7c+8:8^ 

78+0.06 
645+0.06 

6 41+8:8^ 

°-*'-0.07 
6 34+0.06 
°'^'-0.07 
6.20;''"5 

5 97+8:8f 

^■^'-0.06 

6 29+" '*' 
"■^'-0.07 
6 08+" 

5 73+8:8^ 

538+8:8? 



5 58+0.06 

4 76+8:8§ 

-0.07 
r I O+0.06 
= ■"'-0.07 
r oc+0.06 
^'■'"'-0.06 

5 Qfi+0.05 
= ■'"'-0.06 
5.90+"95 



52 03+0.07 
= ■ -0.08 
52 95+0.07 

52 77+0.05 
=^■"-0.10 

52 60+" "* 
■'^- -0.09 
5241+0.07 

57 38+8:8? 

l-r+0.08 

^'■''S-o.io 

52 24+0.08 
-0.10 
52.51+0.08 

52 65+8:8? 

7Z.07_Q iQ 

52 44+0.09 
-0.10 

57 33+0.07 

■ -Q.Q8 

5149+™ 
= ^■^^-0.08 
57 38+0.09 

^^■^''-o.io 
5201+0.08 

51 03+8:8? 

■ -Q.Q8 
50 48+0.08 



4.95+ 



5.86+ 
5.85+ 



49 26+o:o7 

■ -Q.QS 
52 04+''^"' 

51.05+™ 
48.87+g;?;5 

50 54+0.09 
^^■^^-0.10 

f,+0.08 
-0.09 

51.60+g;S 
51-68tS 

51 68+0.08 
51.68„„, 

51 62+""'* 

■ -0.08 



51.36+ 



4.71 
4.91- 
5.03; 
5.70+1 



+0.0 



-0.08 

4 95+0.07 
■ -0.08 

497+0.07 

5 75+8:8i 

-0.08 

5 36+" "7 
■^■^"-0.07 
5 35+0.06 

4 17+8:8? 
*-"-0.07 
c IO+0.07 
'■'"-0.08 
478+0.08 

407+8:85 
3 73+8:8^ 



4 48+u.oa 
^■^°-0.08 

4 16+°"7 
+■'"-0.08 

4 64+""* 
+■ -0.08 

4 63+""* 
+■"^-0.08 
4 95+0.07 
*-'=-0.07 
4 7 1 +0.08 
-0.08 



4 76+"'" 
+■' -0.09 
4 65+0.10 

4 63+8:18 
*-°^-0.10 
4 08+012 
+■ -oil 
4 54+0.10 

4 71+8;?? 

5 47+8:1? 
'■*'-0.11 
4 67+010 
+■"'-010 

4 61+010 

5 00+8:1' 

= ■""-011 
5 37+0.10 
■^■=-010 
531+010 

3 58+8:?? 
=■==-010 
5 00+011 
■'■""-o.io 

4 41+011 

3 58+8:1''' 

3 15+8:15 

^■"-011 
2 94+0.09 

462^8:?? 
+■ -oil 
4.10 

1 

4.39 
3.60 
4.3 1_ 

4 25+o;ii 
+■^■^-011 
4 84+010 
^■''^-OIO 
4 43+011 
+■ -oil 



-0.11 
+0.11 
0.11 



+0.11 



T C7+0. 1 1 

^■'*2-oii 

4 ,6+0.09 

4 03+8:18 
+■"=-011 

4 18+""^ 
+■' -0.06 

3 94+0.09 

434i8;41! 

4 95+8:8? 
+■^=-0.08 
3 04+010 



4 43+0.U8 
+■^=-0.09 

4 10+"'^ 
+■'"-011 



439+0.09 

3 42+8:i8 
=■^^-0.07 
421+010 
+■^'-011 

4 17+0.05 

324+8:8§ 
3 22+8:8^ 



-0.09 
3 74+0.06 
=■ '+-0.07 
, oc+0.06 
'■''=-0.07 
3 91+0.06 
=■"-0.07 
4+0.06 
+-0.07 



3.44+ 



4 00+" "^ 
+■""-0.06 
403+0.05 

4 io+8:8§ 

+■'"-0.09 

4 06+""^ 
+■""-0.06 



I 159+0.024 
-0.027 
1.350;""2'' 

1 304+8:852 

'■="*-0.034 

1 261 +""2* 
'■^"'-0.031 
1 226+""24 

'■^■^"-0.(17^ 
^-0.028 



1.262+ 



1 7O9+0.030 
'■^"^-0.034 
1.254+""'" 
I 334+8:85(i 
'■==+-0.034 

, 7or+0.031 

'■'^"=-0.035 
1.288!";"^^ 

0.965ii 
1 267+0.030 
'■^"'-0.034 
1 167+0.027 

O 867+8:8^8 

746+"^"2* 

"■'+"-0.030 
395+0.025 

"85S 
' -0.034 



1.177+ 



261+""™ 
"■^"'-0.034 
A 07^+0.030 
"■8 '=-0.033 

945+0.026 
"■^+=-0.030 

1 073+0.026 

-0.028 
1 088+0.026 
'■'"'8_o.028 
1 120+0.027 
'■'^"-0.031 
1 O86+0.026 
'■'""'-0.029 



,783+ 



+0.024 
-0.027 
963+0.025 

922+8:858 



842- 



^-0.035 
+0.031 
-0.033 
841+0.025 

860+8:85? 
:8:85? 

' = '-0.035 
8 1 Q+0.030 
""-0.034 
869+2^23'; 
If 



Mil 
-0.034 
+0.025 



^^"^-0^028 
og-1+0.031 
^^^-0.035 
75O+0.029 

:8:83? 
:8:811 

-0.033 
037+0.024 

70 I +8:851 

""-0.036 
533+0.029 
===-0.032 
+0.031 
-0.035 
462+0.032 

545+8:831 
=+=-0.031 

666+""2* 
"""-0.030 
676+0.028 
"'"-0.030 
778+0.028 
'^"-0.032 
679+0.028 
"' -0.031 



9.74+ 



10 11+0.09 
'"■"-0.09 
Q 81+0.03 



9 61+8.85 
^■"'-0.06 
9 79+0.09 
■^-0.09 
9 93+0.09 

9 56+8:8? 
^■="-0.07 
9 74+010 
'■'+-010 
934+0.09 

9 80+8:8? 

^■""-0.06 
9 59+0.09 
-0.09 
936+0.05 

9 65+8:8t 
-S>M 
9 26+0.03 
'■^"-0.03 
9 81+0.09 

9.45+' 



-0.07 
9 ,0+0.05 

9 58+8:8? 

'■="-0.05 
o ^7+0.04 
"■=^-0.04 
8+0.06 
^^-0.06 



+0.05 



9.68+; 
9 24 

'■^+-0.05 

9 26+" "* 
'■^"-0.04 
9 16+0.08 
^■'"-0.07 
9 17+0.05 
'■' '-0.05 



069+""''7 
"■""'-0.090 
nci<; +0.095 



-0 928+" "*3 
"■'■"'-0.076 
-0.356+"™5 



-0 155+"."''7 
"■'==-0.091 
_0.370+"."''7 

_0 2l7+8:8?3 
-0.087 

255+"'"" 
"■^==-0.094 

0.225;""''= 

-0 533+8:852 
"■===-0.082 
01 8+0.097 

-"■""* -0.090 
O 3 1 8+0.086 
-S.Q39 

-0.239 



:8:8i5 

-0.075 
+0.099 

-'^-^^'^ -0.086 

A 3^^+0.090 

A 54Q+0.086 

c+0.090 
■^-0.082 
A no+0.087 

A 9-1 0+0.085 

264 

-0 050+^-^^^ 



0.542 



-0.455+ 



A machine-readable version of Table [A3| is available online (http://www.astro.keele.ac.uk/-dem/clusters.html) or upon request from the first author. 



Table A4. Galactocentric radii and *:-space pa rameters from I Kind 
model fits to 27 GCs from lHarris et al.l j2002l = |Paper]|) 



Name 


Detector 


[kpc] 


Model 


fm,l 


«m,2 


«m,3 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


C002 


HAR02/V 


14.20 


K66 


n n/;Q+0.049 


4 97 1+0120 
+■'"-0.120 


0.416!S;g|7 


C007 

coil 


HAR02/V 
HAR02/V 


9.20 
6.60 


K66 
K66 


0.399+0.048 

354;8:8?S 
"■==+-0.059 


5 372+"^"^7 

5 2i3;8:l?§ 

-0.126 


0.388!»;"|^ 
399+8.0I7 
"■="-0.016 


C017 


HAR02/V 


6.20 


K66 


186+""*- 


5 365+0.064 
= ■="=-0.073 


345+0.005 
"■=^=-0.003 


C021 


HAR02/V 


7.20 


K66 


i98J8:8J1 
"■"°-0.050 


5 106+"'"^ 
= ■'""-0.111 


394+0.017 
"■= -0.015 


C022 


HAR02/V 


5.70 


K66 


0.065+""'"' 
43518:8!? 
"■+==-0.053 


5 5 64+0.067 
= ■="+-0.077 


359+0.011 
"■=='-0.008 


C023 


HAR02/V 


5.80 


K66 


6 3 1 8+0.083 


366+" "'-' 
"■=""-0.009 


C025 


HAR02/V 


8.50 


K66 


1 24+0.055 


5.110+"'5'' 


399+0.017 


C029 
C031 


HAR02/V 
HAR02/V 


21.00 
6.10 


K66 
K66 


0233+8:851 

i8i+8:8S 
0.209+8:855 


5.1 84+8: !;5 

5 694+8:''''^ 
=■"^+-0107 


0393+8:815 

n 774+8:815 
''-™-0.011 


C032 


HAR02/V 


12.50 


K66 


5 323+" 1" 
=■="-0.130 


0433+0.015 


C037 


HAR02/V 


11.90 


K66 


058+8:858 

"■"=°-0.050 


r ^:r7+0.103 

5.652„,i, 


393+8:81? 
"■=^=-0.015 


C040 


HAR02/V 


23.20 


K66 


-0.136+""*5 


4445+0.073 


0366+"."12 


C041 


HAR02/V 


23.20 


K66 


123+8:8^? 

■ -0.058 


5 44o+8:?f 5 
=■++"-0.123 


394+8:819 
-0.015 


C044 


HAR02/V 


20.20 


K66 


-0.077+""*' 


5 354+0.064 


3 52+0.008 

376+8:8'''? 

"■='"-0.012 


ClOO 


HAR02/V 


8.20 


K66 


-0 797+8:8i 
" -"-0.051 


4 230+8:8S 

+■^="-0.101 


ClOl 


HAR02/V 


6.20 


K66 


-0 585+" "*' 

-0 743+8:8fi 
"■'+=-0.050 


4 182+""^* 


0344+0.004 


C102 


HAR02/V 


5.40 


K66 


2.584+8:'^^ 


452+8:8'>'5 

"■+=^-0.014 


C103 


HAR02/V 


5.90 


K66 


-0.036+""" 

-0 426+8:85^ 

-D.94S 


c 7n/;+8:i'l^8 
= ■='"'-0.103 


370+0.013 
"■='"-0.011 


C104 


HAR02/V 


8.90 


K66 


A S7K+0.068 


0357+0.010 


C105 


HAR02/V 


9.10 


K66 


-1 007+""*' 
'■""'-0.052 


7 5 1 7+0.()76 
-0.088 


357+8:88? 
"■=='-0.007 


C106 


HAR02/V 


9.20 


K66 


-0 717+0.046 


c n7fi+0.070 
= ■"=8-0.079 


348+0.00'i 
"■=+"-0.004 


G019 


HAR02/V 


9.00 


K66 


-0 202+8:852 
"■^"''-0.048 


4 467+"^"73 
+■+"'-0.085 


363+0.01 1 
"■="=-0.009 


G221 


HAR02/V 


9.40 


K66 


-0.230+""*' 


5 153+0.063 
= ■'==-0.072 


3 54+0.009 


G277 


HAR02/V 


8.60 


K66 


-0 703+8:8^5 

" -"=-0.043 


5 187+0.060 

5.292+8:85^ 
5 775+8:i')85 
-0.076 


3 50+8:88^ 
"■=="-0.005 


G293 


HAR02/V 


9.50 


K66 


-0.197+""-''" 


379+0.015 
"■="-0.012 


G302 


HAR02/V 


7.90 


K66 


-0 236+8:85' 
"■^="-0.046 


3 55+0.009 
"■===-0.006 


A machine-readable 


version 


of 


Table |A4l 


is available online 



(http://www.astro.k.eele.ac.uk/-dem/clusters.html) or upon request from the first author. 
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